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Abstract! 

The  potential  of  high  intensity  pulses,  in  particular 
socalled  zero-pi  pulses,  is  investigated  for  transmission  of 
submillimeter  waves  through  the  atmosphere.  The  water  vapor 
lines  cause  absorption  which  can  be  partly  reduced  via  satura- 
tion effects. 

In  the  theoretical  part  of  the  paper,  the  inverse  scattering 
method  as  applied  to  selfinduced  transparency  of  a lossless 
medium  is  reconsidered,  physical  interpretations  of  the  mathe- 
matical steps  are  introduced,  and  some  errors  that  appeared  in 
the  literature  are  corrected.  It  is  found  that  the  zero  pi 
pulse  "soliton"  solution  of  a lossless  medium  (no  collisions) 
exists  only  for  pulses  with  the  spectrum  centered  at  the  medium 
line  center,  and  that  pulses  with  phase  reversal,  like  the  zero- 
pi  pulse,  but  with  the  spectrum  put  off  line  center,  break-up. 

A numerical  investigation  of  pulse  propagation  in  the 
presence  of  loss  is  made.  The  analysis  is  carried  out  in  two 
stages.  First,  we  look  at  short  pulse  propagation  with  a carrier 
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frequency  at,  or  near,  the  6.1  cm  1 absorption  line  of  H20 
with  T^  = T2  = 183  psec.  The  line  was  chosen  because  it  lies 
in  the  low  absorption  part  of  the  spectrum  and  permits  propaga- 
tion over  the  order  of  100  km  permitting  a test  of  the  stability 
of  the  computer  program  for  long  distance  propagation.  This  part 
of  the  analysis  shows  that  for  the  carrier  frequency  centered  at 
the  absorption  line  the  loss  experienced  by  a pulse  with  phase 
reversal  may  be  reduced  initially  to  55%  of  the  loss  that 
would  be  experienced  by  a small  signal  pulse  of  the  same  spectrum. 
For  a pulse  spectrum  detuned  so  that  one  of  the  maxima  coincides 
with  the  peak  of  the  absorption  line,  the  reduction  is  greater, 
to  25%  of  the  linear  loss.  As  soon  as  an  appreciable  portion 
of  the  energy  is  depleted,  the  linear  limit  of  absorption  is 
reached.  Typical  intensities  for  this  nonlinear  effect  to  happen 
are  8 MW/cm2  for  a pulsewidth  of  100  psec. 

Contrary  to  some  statements  in  the  literature,  the  zero-pi 
pulse  does  not  possess  a small  signal  limit  because  it  calls  for 
full  inversion  during  part  of  its  evolution.  For  this  reason, 
an  optical  beam  of  finite  cross  section  that  initially  behaves 
like  a zero-pi  pulse  across  its  cross  section  is  subject  to 
diffraction  effects  not  unlike  the  27r-pulse. 

In  the  second  stage  of  the  numerical  analysis  the  two  H20 
lines  adjacent  to  the  methyl  fluoride  (CH3F)  laser  line  (496  u ) 
and  D20  laser  line  (385  \i)  are  investigated.  The  level 
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degeneracy  is  taken  into  account.  The  dipole  moments  are  of 
the  same  order  as  that  of  the  6.1  cm  ^ line  and  therefore 
the  intensity  levels  are  comparable.  The  degeneracy  counteracts 
somewhat  the  nonlinear  loss  reduction,  but  not  greatly. 

Linear  propagation  of  the  CH3F  line  and  D20  line 
achieves  greater  propagation  distances  (0.57  km  and  1.13  km 
respectively  for  the  1/e  decay)  than  nonlinear  propagation  at, 
or  near,  the  respective  absorption  lines. 

The  report  contains  many  computer  plots  which  aid  in  the 
physical  understanding  of  nonlinear  propagation. 


4 


Short  Pulse  Propagation  in  the  Submillimeter  Region 

Introduction 

The  use  of  laser  radar  as  opposed  to  microwave  radar  in 
tactical  weaponry  has  the  advantage  of  greater  resolution.  Ten 
micrometer  CO 2 laser  systems  are  under  development  for  radar 
applications.  This  wavelength  does  not  lend  itself  to  all  weather 
operation  and  can  be  countermeasured  by  dispersing  aerosols.  Sub- 
millimeter  waves  would  not  be  so  susceptible  to  fogs  and  aerosols 
but  would  still  afford  higher  resolution  than  microwaves.  In  the 
submillimeter  region  molecular  absorption  becomes  a severe  problem, 
particularly  absorption  by  pure  rotational  transitions  in  water. 
Fig.  1 gives  a sea  level  1 km  transmission  spectrum  in  this 
wavelength  region. 

It  is  known  that  intense  pulses  much  shorter  than  T2  may 
propagate  through  a resonantly  absorbing  medium  without  experi- 
encing loss.  This  phenomenon  is  called  self -induced  transparency. 

A pulse  with  an  electric  field  of  "area" 


}co 

E dt  = e 

— OO 


2tt 


(1) 


where  y is  the  dipole  moment  and  ft  is  Planck's  constant, 
can  propagate  through  a inhomogeneously  broadened  medium  without 
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attenuation.  The  pulse  is  so  intense  that  it  inverts  the  two 
level  system,  initially  in  the  ground  state,  and  then  returns  it 
to  the  ground  state  so  that  no  net  absorption  occurs. 

If  one  intends  to  utilize  this  phenomenon  for  long  distance 
propagation  of  intense  pulses,  one  quickly  finds  that  the  natural 
spreading  of  the  beam  diameter  via  diffraction  changes  the  "area" 
of  the  electric  field  assigned  to  different  points  in  the  beam 
cross  section  from  its  value  of  2n.  When  the  beam  spreads  the 
decrease  in  intensity  causes  a decrease  of  the  field  amplitude, 
and  hence  a decrease  of  the  "area"  to  less  than  2tt.  When  this 
happens,  an  appreciable  portion  of  the  population  of  the  two 
level  system  is  left  in  the  upper  level  thus  causing  loss  and 
a rapid  absorption  of  the  optical  pulse. 

Socalled  zero-ir  pulses  (Lamb,  1974)  have  a phase  reversal 
so  that  0 of  Eq.  (1)  integrates  to  zero.  Offhand  one  might 
expect  that  zero-ir  pulses  would  be  affected  less  by  diffraction 
effects  than  2tt  pulses.  The  question  then  arises  whether  the 
potential  for  low-loss  propagation  of  zero-ir  pulses  may  be 
utilized  for  transmission  of  mil lime ter -wave  pulses  through  the 
atmosphere,  decreasing  the  effect  of  water  vapor  absorption. 

Also,  computer  solutions  have  shown  that  zero-pi  pulses 
evaluated  for  nondegenerate  two  level  systems  propagate  with 
little  loss  in  a level  degenerate  system  (Lamb,  1974).  There- 
fore, particular  attention  has  been  devoted  to  the  study  of 
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propagation  of  the  zero-pi  pulse. 

To  date,  only  the  zero-pi  pulse  with  carrier  frequency 
tuned  to  the  center  of  a symmetric  line  of  an  inhomogeneous ly 
broadened  system  has  been  treated  in  the  literature.  Let  us 
briefly  review  the  principal  findings  and  explain  their  importance 
for  the  study  of  low  loss  propagation. 

The  zero-pi  pulse  changes  phase  by  180°  in  its  "life-time". 
The  individual  positive  areas  and  negative  areas  of  the  electric 
field  as  defined  by  the  integral  over  time  (1)  are  of  the  order 
of  2tt  however,  i.e.  the  E-field  never  becomes  small  in  the 
strict  sense  of  the  word.  Contrary  to  popular  belief  (Lamb,  1974) 
a small  signal  analysis  is  never  applicable  because  the  medium  ex- 
periences full  inversion  and  returns  to  the  ground  state  within  the 
lifetime  of  a pulse.  (Remember,  this  is  a discussion  for  the  colli- 
sion free  case.  The  field  has  a longer  time  to  act  if  the  amplitude 
is  reduced.)  Therefore,  a zero-pi  pulse  of  a given  duration  ceases 
to  be  a zero-pi  pulse  when  reduced  in  amplitude,  i.e.  it  will 
not  leave  the  medium  in  the  ground  state,  just  like  a pulse 
scaled  down  in  amplitude  from  a 2n  pulse.  Yet,  the  situation 
is  not  as  bad  for  the  zero-pi  pulse.  Because  of  its  phase  re- 
versal, part  of  the  "damage"  done  in  one  part  of  the  pulse  is 
undone  during  the  phase  reversal  part,  less  population  is  left 
in  the  upper  state  after  passage  of  a "reduced  zero-pi  pulse" 
than  after  the  passage  of  a "reduced  2n  pulse".  One  confirma- 
tion of  this  fact  are  computations  (Lamb,  1974)  on  a system  with 


a degenerate  level  in  which  the  differences  of  dipole  moment 
cause  different  time  evolutions  of  the  individual  systems.  A 
degenerate  system  causes  less  loss  to  a zero-pi  pulse  than  to  a 
2tt  pulse.  For  this  reason,  the  effect  of  diffraction  is  to  be 
expected  to  be  less  on  the  zero-pi  pulse. 

The  present  work  is  an  analysis  of  propagation  of  the  zero-pi 
pulse,  with  carrier  frequency  on  line  center  and  off  line  center, 
in  the  presence  of  collisions.  It  consists  of  an  analytic  part 
and  a computational  part.  The  analytic  part  is  concerned  with  an 
application  of  the  Inverse  Scattering  Method  (ISM)  to  the  evalu- 
ation of  the  zero-pi  pulse.  Previous  work  in  the  literature 
(Lamb,  1973,  1974)  has  not  made  mention  of  a zero-pi  pulse 
solution  in  the  offresonant  case  (carrier  frequency  detuned 
from  the  center  of  a symmetric  line).  Offhand,  one  would  expect 
that  in  a medium  with  multiple  lines  lower  loss  is  achieved  when 
the  carrier  frequency  of  the  pulse  is  placed  between  the  lines. 
Therefore,  it  was  necessary  to  apply  the  analytic  theory  of 
zero-pi  pulses  to  this  case. 

A reader  of  the  literature  on  ISM,  who  is  not  a professional 
mathematician,  tends  to  be  dazed  by  the  seeming  complexity  of  the 
subject.  This  was  the  case  with  the  author  of  this  report.  We 
perceived  it,  therefore,  to  be  our  first  task  to  cast  the  termin- 
ology developed  in  the  ISM  into  language  familiar  to  the  physicist 
and  electrical  engineer.  This  topic  is  covered  in  detail  in 
Appendix  I.  Suffice  it  to  state  here  that  we  succeeded  in  sim- 
plifying greatly  the  derivation  of  the  zero-pi  pulse. 
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The  computational  part  studies  propagation  of  pulses  of 
zero  area  in  the  presence  of  collisional  deexcitation.  The  loss 
of  a zero-pi  pulse  is  compared  with  the  loss  that  would  be  ex- 
perienced by  an  equivalent  low  intensity  pulse.  By  this  we  mean, 
a pulse  with  the  same  spectrum  and  frequency  dependent  loss 
matching  the  linear  transmission  characteristic  of  the  medium. 

The  investigation  is  divided  into  two  parts.  First  we  ignore 
the  orientational  degeneracy  and  study  nonlinear  propagation  in 
a medium  with  loss.  The  parameters  of  the  medium  are  picked  so 
as  to  model  the  water  absorption  line  near  6.1  cm  ^ where  the 
loss  is  relatively  small  and  long  distance  propagation  can 
occur.  This  is  one  test  on  the  stability  of  the  computer  pro- 
gram. This  part  of  the  study  gave  a great  deal  of  insight  into 
the  possibility  of  loss  reduction  by  nonlinear  effects.  It  was 
found  that  the  loss  of  a pulse  of  100  psec  duration  of  8MW/cm 
peak  power  could  be  reduced  to  about  20%  of  the  loss  experienced 
by  a linear  medium.  Of  course  this  reduction  lasts  only  for  a 
few  1/e  attenuation  lengths  after  which  the  loss  is  that  of 
a linear  medium.  It  was  also  found  "empirically"  in  this  in- 
vestigation that  zero-pi  like  pulses  tuned  off-resonance  tend 
to  break  up  under  lossless  propagation.  This  prompted  the 
analytic  study  to  look  for  a general  proof  that  detuned  zero-pi 
pulses  do  not  exist  which  is  reproduced  in  the  appendix. 

The  second  part  of  the  study  is  concerned  with  a realistic 
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modeling  of  the  propagation  of  pulses  at  resonance  with  the 
absorption  lines  adjacent  to  the  495  y line  of  the  CH3F 
and  the  385  y line  of  the  DjO  laser.  Here  we  took  the 
orientational  degeneracy  into  account  finding  that  it  changes 
the  propagation  characteristics  only  very  little.  The  penetra- 
tion depths  at  these  frequencies  are  discouragingly  low  however. 
Greater  propagation  distances  are  achieved  by  small  signal  pro- 
pagation of  the  CH^F  and  D20  lines  proper,  because  they  lie 
relatively  close  to  local  absorption  minima. 


A large  amount  of  literature  has  evolved  in  recent  years 
(see  References)  on  the  Inverse  Scattering  Method  (ISM)  which 
yields  exact  solutions  to  a class  of  nonlinear  partial  differ- 
ential equations  in  two  dimensions  (time  and  one  space  variable). 
The  equations  of  self-induced  transparency  (S.I.T.)  are  a member 
of  this  class,  and  the  zero-pi  pulse  is  a solution  of  these 
equations  that  has  been  derived  by  means  of  the  method  (Lamb,  1974) 
It  seemed  appropriate,  therefore,  to  devote  some  effort  to 
the  exploration  of  the  method.  It  became  soon  apparent  that  the 
research  on  the  I.S.M.  had  been  conducted  exclusively  by  mathe- 
maticians and  had  not  been  assimilated  by  physicists  and  engineers, 
like  circuit  theory  had  been  adapted  in  the  thirties  after  its 
original  development  by  mathematicians.  For  this  reason  we  set 
out  to  develop  an  understanding  of  the  I.S.M.  Appendix  I is  the 
outcome  of  this  study.  Here  we  summarize  the  principal  issues. 

The  flow  chart  of  Fig.  2 will  help  the  reader  in  following  the 
development. 
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The  I.S.M.  associates  with  the  original  nonlinear  differen- 

/ 

tial  equation  a linear  (quantum  mechanical)  scattering  problem, 
the  solution  of  which  provides  the  general  solution  of  the  non- 
linear differential  equation.  This  association  is  done  ad  hoc, 
no  general  methods  exist  to  find  the  linear  scattering  problem. 

In  the  S.I.T.  problem  Lamb  (1973)  found  the  linear  scattering 
equations,  the  socalled  Zakharov  Shabat  equations,  after  a series 
of  ingenious  transformations  of  the  Bloch  equations.  Because 
S.I.T.  is  a quantum  mechanical  problem  (so  far  as  the  medium 
description  is  concerned)  one  would  expect  that  the  linear 
scattering  problem  is  naturally  imbedded  in  the  inverse  scat- 
tering method  applied  to  S.I.T.,  in  other  words,  is  part  of  the 
defining  equations.  This  has  not  been  generally  recognized  and 
only  McLaughlin  and  Corones  (1974)  have  made  the  connection,  but 
not  explicitly.  We  showed  that  the  linear  scattering  equations 
of  S.I.T.  are  the  equations  of  the  two-level  system  before  they 
have  been  cast  into  the  Bloch  equation  form. 

Next  we  study  the  solutions  of  the  Zakharov-Shabat  equations. 
Here  we  establish  the  analogy  of  these  equations  with  those  of 
the  parametric  oscillator—  a well  known  problem  in  Optical 
Electronics  (Yariv) . This  analogy  is  helpful  in  locating  the 
poles  in  the  complex  plane  of  the  eigenfunctions  of  the  Zakharov- 
Shabat  equations  that  give  the  "soliton"  solutions  of  S.I.T. 

These  poles  must  be  invariants  of  the  S.I.T.  equations  if  the 
inverse  scattering  method  is  to  be  applicable.  We  give  a proof 
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of  the  invariance  following  a method  of  Ablowitz  et  al.  (1974a) . 
The  invariance  leads  to  a set  of  equations  that  are  shown  to  be 
identical  with  Hilbert  Transforms  of  the  Bloch  Equations.  The 
invariance  may  be  utilized  to  obtain  the  spatial  evolution  of 
the  pulse  determined  at  one  value  of  the  spatial  coordinate; 
the  entire  spatial  dependence  is  contained  in  the  residues  of 
the  poles,  the  location  of  the  poles  is  invariant.  Because  the 
residues  may  be  shown  to  be  independent  of  the  time  variable  t 
evaluation  may  occur  in  the  limit  t -»  which  leads  to  great 
simplification. 

In  this  manner  we  obtain  the  2 n -pulse  and  zero-pi  pulse. 

We  can  also  show  that  no  soliton  solution  with  zero  area  exists 
off  resonance  (carrier  not  centered  with  absorption  line).  In 
this  connection  one  should  emphasize  the  shape  of  the  spectrum 
of  a zero-pi  pulse  (or  a zero  degree,  the  terminology  used  by 
Grieneisen  et  al.,  1972,  1973,  who  considered  small  signal  pro- 
perties of  pulses  with  phase  reversal) . The  spectrum  of  a zero- 
pi  pulse  has  a zero  at  the  carrier  frequency.  This  accounts  for 
the  low  loss  observed  with  zero-degree  pulses  (Grieneisen  et  al., 
1972)  when  propagating  through  a medium  whose  absorption  line 
was  centered  with  the  carrier  frequency. 
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II.  The  Zero-Pi  Pulse 

The  zero-pi  pulse  with  a carrier  frequency  centered  with 
the  symmetric  line  of  a two  level  system  distribution  has  been 
derived  by  Lamb  (1974)  from  the  Inverse  Scattering  Method.  It 
is  rederived  in  the  appendix  and  has  the  following  form  in  the 
notation  of  the  appendix: 


f - 8S 


cosh  28(x  - xQ)cos  2a(x  - x^  - — sinh  2B(x  - xQ)sin  2a(x  - t -^ ) 

a 

cosh2  28 (t  - t ) + — sin2  2a (t  - x.  ) 

0 9 X 

a 


(2.1) 


a and  8 are  adjustable  constants.  F is  the  normalized  E 
field 
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where 


oj  - carrier  frequency 
N - particle  density 
■fi  - Planck's  constant 

e - dielectric  constant  of  (nonresonant  part  of)  medium 


Further,  x is  the  normalized  space-time  variable  which  refers 
to  a frame  translated  along  the  x-axis  at  the  speed  of  light 


x = 


(2.4) 


The  normalized  spatial  variable  is 


z = ftx/c 


in  terms  of  which  one  may  write 


= — z <■ 


4 (?  - a)2  + e2 


ax,  = — z < — ~ > 


4 (C  - a)2  + e2 


(2.5) 


(2.6) 


(2.7) 


The  pointed  brnckrts  i nd  i on  to  .in  nvt'rngo  over  nil  two  lovol  systems 


. 
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of  different  transition  frequency  (inhomogeneous  broadening). 

In  the  limit  of  no  inhomogeneous  broadening,  5 can  be  set 
equal  to  zero  in  the  above  equations  and  the  pointed  brackets 
can  be  omitted.  This  is  the  limit  of  particular  interest, 
because  the  lines  of  water  vapor  in  the  centimeter  and  milli- 
meter region  in  the  atmosphere  are  homogeneously  broadened. 

In  this  limit  one  finds  that  the  zero-pi  pulse  consists 
of  an  envelope  traveling  at  less  than  the  speed  of  light,  and 
a substructure  with  a phase  velocity  greater  than  the  speed  of 
light.  The  inverse  speeds  are  disposed  symmetrically  around 
1/c. 

[ ( Some  features  of  the  zero-pi  pulse  are  apparent  from  (2.1) 

and  also  follow  from  its  derivation  via  the  Inverse  Scattering 
Method.  The  zero-pi  pulse  may  be  considered  to  be  constructed 
by  a limit  of  two  2tt-pulses  that  propagate  at  the  same  speed, 
are  in  antiphase  and  are  made  to  overlap.  This,  of  course,  is 
a construction  in  the  "indirect  space"  of  the  linear  scattering 
protlem.  not  in  the  "direct"  space  of  the  nonlinear  differential 
equation,  in  which  superposition  does  not  hold.  An  indication 
of  this  "superposition"  is  the  fact  that  the  sech2  of  the 
two-pi  pulse  has  the  same  scale  parameter  (5,  which  also 
appears  in  the  amplitude  of  the  pulse  in  front  of  expression 
(2.1).  This  means  that  the  individual  "areas"  of  the  zero-pi 
pulse  are  of  the  order  of  2 tt  and  not  small.  Therefore,  the 


r 
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zero-pi  pulse  shares  with  the  2n  pulse  the  property  that, 
strictly,  a small  signal  limit  does  not  exist.  In  both  cases, 
appreciable  inversion  occurs  during  the  time  evolution  of  the 
pulse.  Therefore,  the  zero-pi  pulse  is  susceptible  to  diffrac- 
tion effects  just  like  the  2 tt  pulse.  In  other  words,  spreading 
via  diffraction  that  changes  the  field  intensity  while  not 
changing  the  time  evolution  will  change  the  positive  and  nega- 
tive portions  of  the  pulse  in  different  parts  of  the  beam  and 
make  it  become  different  from  the  ideal  zero-pi  pulse.  It  is 
important  therefore  to  study  propagation  of  intense  pulses 
other  than  zero-pi  pulses  by  computer  in  order  to  ascertain  the 
potential  of  intense  pulse  propagation  through  the  atmosphere. 

In  the  appendix,  an  investigation  is  made  of  the  existence 
of  zero-pi  pulses  with  a carrier  detuned  from  line  center.  The 
requirement  that  the  contributions  of  the  two  poles  propagate 
at  the  same  speed  imposes  a constraint  which  is  shown  to  be 
equivalent  to  constraining  the  spectrum  to  be  centered  with  the 
material  line. 
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Part  II 

Computation 

III.  The  Starting  Equations  for  Computation 

For  the  computation  of  pulse  propagation  with  loss,  it 
is  necessary  to  use  the  Bloch  equations,  or  density  matrix 
equations.  The  loss  is  introduced  phenomenologically  in  terms 
of  the  "longitudinal"  and  "transverse"  relaxation  times  T^ 
and  T in  the  slowly  varying  envelope  approximation  of  an 
electric  field,  E(x,  t)  exp  - i (<ot  - kx)  , these  are  (Kryukov, 
Letoknov,  1970)  (compare  (1.10),  (1.15)  and  (1.16)  of  the  ap- 
pendix) 


' 

_3_  + 1 _3_ 
,3x  c 3t. 


E (x , t) 


i2ir  — P(x,  t) 
c 


(3.1) 


— + — i Aw 

at  T_ 


P ( X , t)  = -1^-  ■ N (x,  t)  E(x,  t) 


(3.2) 


a N(x,  t)  - N . 

-i-  N (x , t)  + — = 1 

3 1 t.  -n 


E*  (x,  t)  P(x,  t)  - E (x,  t)  P*  (x , t) 


(3.3) 


with 


All)  = W - U) 

o 


(3.4) 
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where  u>o  is  the  resonance  frequency  of  the  two-level  system. 

We  use  here  cgs  units.  N(x,  t)  is  the  difference  of  the 
population  densities  in  the  lower  and  upper  levels,  N(x,  t)  = 
N(pll  " p22}  of  the  aPPendix;  P(x,  t)  is  the  polarization 
density,  p (=  p12  of  the  appendix)  is  the  matrix  element, 

T2  is  the  dephasing  time,  is  the  energy  changing  collision 

time;  Nq  is  the  equilibrium  population  density  difference. 

In  the  case  of  more  than  one  line  we  have 


_9_  + 1 J_ 
ax  c at 


E(x,  t) 


E i2ir  - P.  (x,  t) 
j c 3 


(3.5) 


iAa) . 
3 


Pj(X,  t) 


II  (x,  t)  E (x , t) 


(3.6) 


— N, (x, 
3 1 3 


t)  + 


N.(x,  t)  - 


i 

h 


E*(x,  t)  Pj(x,  t)  - E (x , t)  P * ( x , t) 

2 


(3.7) 


In  the  case  of  one  single  line,  we  normalize  as  follows: 


t 

in 

units 

of 

fco 

X 

in 

units 

of 

ct 

( 

E 

in 

units 

of 

Eo 

P 

in 

units 

of 

Po 
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N in  units  of  N 


where 


tQ  is  arbitrary 


With  this  normalization,  (3.1) -(3.3)  become 


_3_ 

3x 


+ 


E(X, 


3,1 

— + lAo) 

3 1 T2 


— N(x,  t)  + N-(-X->.  — } 
3t  Tx 


t)  = - <xP(x,  t) 


P(x,  t)  - - N(x,  t)  E (x,  t) 


~ 1 = E* (x , t)  P(x,  t)  + E (x 

2 


(3.8) 


(3.9) 


t)  P*(x,  t) 

(3.10) 


where 


20 


2itw  t 2N  I y I 2 
o o'1 


(3.11) 


In  the  case  of  many  lines,  the  normalization  is  slightly  different. 
If  we  define 


Po  “ = l*j 


(3.13) 


and  express  y^  in  units  of  yQ,  then  we  recover  the  following 
normalized  equations 


— + — E(x,  t)  = 
3x  at 


- E a.  P.(x,  t) 

j 3 3 


(3.14) 


— + — iAu) . P.  (x,  t)  = - y.N.  (x,  t)  E (x,  t)  (3.15) 

at  t2_.  3J  3 33 


N.(x,  t)  + 


N.  (x,  t)  - 1 


Tlj  " ^ 


E*(x,  t)  Pj  (x,  t)  + E(x,  t)  Pj*(x,  t) 


(3.16) 


where 


a j = 


2Tra)t  2y  y.  N . 
o o ] 03 


(3.17) 


and 


P^(x,  t)  in  units  of  PQ  = V*jN0j 
Nj(x,  t)  in  units  of  N0 j • 
Numerical  values  are 


Eo  I = — " 6.5821  x 10-8  (touo)_1  ^£5. 

u t cm 

o o 


2irut  2N  |u  |2  4 t 2N  y 

a = P-  = 2.59  x 104  SL  -°-  ° - 

* X 


I s I = — E X H = 2.65  x 10-3  E 2 Watts/cm2 

4 TT 


where 


uo 

is 

in 

e8 

(leA  = 4.80  Debye) 

fco 

is 

in 

sec 

N 

o 

in 

-3 

is 

cm 

X is  in  cm 


(3.18) 

(3.19) 

(3.20) 


is  in  Volts/cm 
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IV.  Differencing  the  S.I.T.  Equations 

In  this  section  we  develop  the  differencing  scheme  that 
has  been  used  in  the  computer  program  to  evaluate  the  pulse 
propagation.  The  S.I.T.  equations  in  normalized  units  are  given 
by  (3.14)  through  (3.16).  If  the  polarization  is  now  split  into 
real  and  imaginary  parts, 


Pj(x,  t)  = PRj  (x,  t)  + iP^U,  t). 


(4.1) 


then  (3.15)  and  (3.16),  separated  into  real  and  imaginary  parts, 
can  be  written  in  matrix  form 


N . 

3 

P . 

3t 

R3 

PIj. 

Tu 

Vr 


-Vi 


U , E_ 
3 R 


2 j 

A co 


Vi 


-Aco 

1 

T 

2j 


N . 

3 

1 

Tii 

PRj 

+ 

0 

pIj. 

0 

.2) 


If  we  define 


N . 
3 


Rj 

> 

Ij 


(4.3) 


^ i 


23 


A. 

3 


1 

T 

UjER 

“jEi 

ij 

-“jER 

1 

T 

2j 

-Au) 

-UjEj 

Aoo 

1 

•p 

lljJ 


(4.4) 


b.  = 
D 


VTlj 

o 

0 


(4.5) 


Then 


— y = A • y + b 
3 1 


(4.6) 


where  the  j subscripts  are  omitted. 

The  strategy  is  to  solve  the  polarization  and  population 
equations  at  a given  point  in  space  for  all  time,  and  then  in- 
crement the  space  point  and  again  solve  for  the  populations  in 
time.  The  time  difference  equations  are 


- 


(4.7) 


(!  - T *i) 


-1 


y_  + — A • y + A tS 
■‘o  2 ° ° 
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?n  + 1 - (!  - 7 \ + 1. 


-1 


- 1 * 7 lSn  - 1 • f„  - 1 


+ 2Sn  ’ ?n> 


+ 2Atb] 


(n  > 1) 


where  yn  and  denote  the  values  at  the  n-th  time  step. 

The  electric  field  is  solved  for  in  the  moving  frame.  If 
we  define 


t = t - x 


and 


z = x 


(4.8) 


as  the  new  variables,  then  (compare  Section  1 of  Appendix)  the 
derivative  3/9 t becomes  3/9 x and  the  "convective"  derivative 
of  E becomes  3/3z  at  constant  x. 


— E(z,  t)  = E a.  P. (z,  t) 


3z 


3 3 


(4.9) 


This  equation  can  be  solved  numerically  using  a predictor- 
corrector  scheme,  where  the  predictor  is 


E1  + 1 = E1  + An  I a -P . 1 

j 3 3 


(4.10) 


and  the  corrector  is 


^ + 1 = — .i,  J 

2 j -3-i 


E*  ' * = E*  - " + ^ [Z  a4(?i  " 1 + 2P . 1 + P1  + 1)] 


(4.11; 


The  algorithm  works  as  follows:  first  the  field  E1  + } is 

n + l 

predicted  according  to 


E 


i + 1 
n + 1 


3n  + 1 + VVn  + 1> 


(4.12) 


The  polarization  and  population  is  then  computed  by 


+ 

+ 


1 

1 


At  51  + 
— A 
2 n + 


+ 1 
- 1 


+ 


At 

2 


+ 1 
- 1 


+ 1 
- 1 


+ 2 


l + 

n 


% + ')  + 24tBo)  • 


(4.13) 


This  is  done  for  all  the  transitions  yielding  values  for  the 
polarization.  With  the  new  polarization  values,  the  field  can 
be  updated  according  to 


„i  + 1 

Jn  + 1 


= E 


,i  - 1 


‘n  + 1 


2 


(p  . ) ~ ^ 

' j'n  + 1 


2(P.)1 
3 n 


+ 1 


+ (P^ 


Vn  + l)  ) * 


(4.14) 


From  the  corrected  value  of  the  field,  the  populations  and 


polarizations  are  recomputed 
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The  code  was  originally  intended  to  use  a quadrature  of  the 

form 


rn  + 1 


= y, 


n 


^ (A 


- 1 * - 1 + 4 


n 


‘ n 


+ A 


n + 1 


rn  + 1J 


However,  this  scheme  when  applied  to  the  field  and  populations 
was  found  to  be  unstable. 


V.  Absorption  Coefficients 

In  this  section  we  give  the  definition  of  the  absorption 
coefficient  and  relations  which  will  be  relevant  in  the  next 
sections.  The  normalized  wave  equation  is  (3,14) 

E(x,  t)  = - z a P (x,  t).  (3 

j J J 


± + -L 

9x  3t 


The  normalized  absorption  nonlinear  coefficient  k is  the  loss 
of  energy  per  unit  normalized  length 


t oo  

J I E | 2dt  f “ | E | 2dt 

“ 00  J _oo 


2dt 


•»  i: 


(E^Pj  + EP*)dt 


K 


(5.1) 
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One  can  define  a linear  absorption  coefficient  by  solving  for 
the  linear  polarization  assuming  that  the  normalized  population 
keeps  its  small  signal  value,  and  then  computing  ic  according 
to  (5.1).  The  ratio  of  the  nonlinear  absorption  coefficient 
to  the  linear  absorption  coefficient  is  a measure  of  the  non- 
linearity of  the  medium  response.  In  order  to  get  from  normalized 
to  absolute  absorption  coefficients,  one  divides  by  (ctQ) . 


VI.  Fourier  Transform 

The  Fourier  transform  of  a function  f(t)  is 
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VII.  Results  - One  Line,  Resonant  Case 

We  shall  attack  the  numerical  solution  in  stages.  We  study 
first  2 7T  pulse  propagation  with  no  loss  and  check  the  computer 
solution  against  the  analytic  solution.  In  this  part  of  the 
problem,  the  choice  of  numerical  values  for  the  parameters  is 
not  crucial.  Indeed,  the  matrix  element  y and  the  number 
density  Nq  jointly  establish  distance  and  time  scales.  A 
change  in  the  values  of  y and  Nq  is  accomodated  fully  by  a 
change  of  the  numerical  values  of  distance  and  time  on  a com- 
puter printout  made  for  a particular  initial  choice  of  y and 

N . 
o 

When  collisions  are  introduced,  they  introduce  a rate 
constant  that  must  be  normalized.  A change  of  y and/or  Nq 
calls  for  a change  of  normalized  rate  constant,  a new  independent 
parameter  has  been  introduced. 

The  same  statements  hold  for  the  lossless  0-tt  pulse  and 
0-tt  pulse  propagation  in  the  presence  of  collisions  (loss)  . 

Hence  one  must  be  careful  to  use  appropriate  numerical  values 
for  the  physical  problem  to  be  solved.  All  this  is  rather  ob- 
vious, but  there  is  additional  complication.  H20  is  an  asym- 
metric rotor  molecule  and  all  levels  of  interest  are  orienta- 
tionally  degenerate;  i.e.  for  a given  angular  momentum  J, 
there  are  2J  + 1 possible  M eigenvalues  ( ’’pro jection  of 


J 
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angular  momentum"  into  the  spatial  axis  around  which  the  angular 
momentum  has  been  quantized) . Each  value  of  M gives  a dif- 
ferent dipole  moment.  Thus,  a given  pair  of  energy  levels 
interacting  nonlinearly  with  an  electric  field  must  be  treated 
as  2J  + 1 level  pairs  where  J is  the  quantum  number  of  the 
lower  level.  Each  level  interacts  with  the  E-field  with  its 
own  value  of  y and  thus  ideal  2ir-pulses  or  0— it  pulses  do 
not  exist. 

In  order  to  break  down  the  problem  into  stages  of  increasing 
difficulty,  we  consider  first  the  case  of  no  M-degeneracy . We 
shall  first  pick  values  for  the  line  strengths  and  widths  from 
the  computer  print -out  of  S.  A.  Clough  based  on  work  by  L.  S 


Rothman  and  R.  A.  McClathey 


i,  15,  2616,  1976) 


in  the  low  absorption  regime  around  6 cm”  to  test  how  the 
program  would  handle  long  distance  propagation.  These  compu- 
tations will  be  useful  in  isolating  the  problems  inherent  in 
2— tt  and  0-tt  pulse  propagation  in  the  presence  of  loss  and 
also  ascertain  the  capability  of  the  program  to  handle  long 
distance  propagation.  Then  we  shall  pick  the  more  realistic 
cases  of  nonlinear  propagation  near  the  two  submillimeter  lasers, 
methyl  fluoride  at  496  y and  D20  at  385  y,  using  H20 
line  parameters  for  the  two  absorption  lines  nearest  the  per- 
tinent laser  line.  Here,  the  full  M-degeneracy  of  the  lines 
will  be  taken  into  account. 
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Table  7.1  gives  the  line  parameters  used  for  the  6.1  cm-1 
propagation. 

Table  7.1  Data  for  6.1  cm-1  line  of  H20. 
frequency  u>  = 1.2  * 1012  sec'1 

relaxation  times  ^ = t2  = 0.094  cm-1  or  185  psec 

matrix  element  u = 0.2204  debye  = 4.58  x io'2  e8 

energy  of  lower  level  131.2  cm”1 

quantum  numbers  J,  K K+1 

lower  level  313 
upper  level  220 

linear  absorption  1.135  x io”6  cm”1 

number  density  in  lower  level  2.77  x io11 

a = 9.18  x io”7 

Eq  = 1.44  x io^  V/cm 

PQ  = 5.50  x watts/cm2 
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Run  1:  Lossless  2 it  pulse  propagation. 

In  the  first  run  we  consider  propagation  of  a 2tt  pulse 
on  resonance,  comparing  the  analytical  results  with  the  numerical 
results  as  a check  on  the  computer  program.  In  this  problem  we 
take  the  following  parameter  values: 


tQ  = 100  psec  A-1  = 6.1  cm"1 

N = 2.77  x 1011  cm-3  Ao)  = 0 
o 

U = 4.58  x 10-2  eft  Tx  = T2  = 1050  sec 

and  for  time  and  space  steps  we  have 


min 

0 

3 x 

lO^cm 

fcmin 

t 

-1  nsec 

1 nsec 

max 

Ax  = 

3 x 

103  cm 

max 

At  = 

6.67  psec 

nx  = 

100 

steps 

nt  = 

300  steps 

(7.1) 


(7.2) 


The  input  field  is  2 sech  t in  normalized  units.  The  time 
axis  was  shifted  on  each  space  step  such  that  the  average 

Z ti((Re  E)?  + (Im  E) ? ) At 

Z ((Re  E)|  + (Im  E) ? ) At 
i 1 1 


<t> 


(7.3) 
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was  kept  between  the  centermost  two  time  bins.  The  column 
marked  R1  in  Table  7.2  gives  the  distance  in  cm  at  which 
the  remaining  quantities  defined  below  are  evaluated. 

AREAR  = E (Re  E) . At  (7.4) 

i 1 

ENERGY  = E (E*E) ^ At  (7.3) 

VELOCITYC  = 1 - 0<t>/3x)_1  (7.6) 

E a.  I ( E*P j + EP^*^  At 

KAPPAl  = 2 i (7.9) 

E (E*E)i  At 

The  summation  over  j in  KAPPAl  is  over  the  different  transi- 
tions, in  this  case  there  is  only  one  considered.  AREAR  is  the 
trapezoidal  rule  integration  for  the  real  area  of  the  electric 
field  and  should  be  equal  to  2ir  for  the  entire  problem.  The 
error  is  seen  to  be  of  order  0.1%  or  less.  The  pulse  energy 
rises  slowly,  with  an  error  of  0.28%.  The  pulse  velocity  is 
the  speed  of  light  to  within  about  one  part  in  106,  the  com- 
puted value  of  VELOCITYC  is  accurate  to  about  0.7%,  rendering 
the  absolute  velocity  accurate  to  7 parts  in  109.  The  ab- 
sorption coefficient  analytically  is  exactly  zero,  the  error  is 
such  that  the  growth  in  energy  is  of  the  order  of  one  part  in 


33 


3®IS  PASS  IS  BEST  QUALITY  FBACIICABLi 
TOO*  COPY  PURKISHEE  TO  DOG 


Table  7.2 


Run  1.  The  2it -pulse,  error  check. 
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10  over  3 x 10  cm  propagation. 

Basically , the  numerical  errors  introduced  are  of  order 

-3  -2 

10  to  10  , which  are  acceptably  small  for  the  calculation. ' 

That  is  to  say,  the  numerical  accuracy  is  of  order  99.0%  to 
99.9%. 
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Run  2:  Lossless  Ott  Pulse  Propagation 


In  this  run  we  use  the  same  parameters  as  given  in  (6.1) 
and  (6.2),  however,  the  input  field  is  now 

t. 


cos (t/t, ) cosh  (t/t  ) - — sin  (t/t,)  sinh(t/t  ) 
i g , r 0 


E (0 , t)  = — 
fce 


cosh2  — + (t,/t  )2  sin2  (t/t,)  sinh2  (t/t  ) 

^ 1 C I.  0 

e 


in  normalized  units.  Vaxues  of  t and  t^  are 

e f 


(7.8) 


t0  = 10  10  sec 


t^  = 10~10  sec. 


(7.9) 


In  Table  7.3  we  give  the  same  values  of  the  parameters  as  listed 
in  Table  7.2,  now  for  the  case  of  0-tt  pulse.  The  pulse  energy 
is  within  0.2%  of  its  analytical  value.  The  velocity  is  correct 
to  within  1%  in  terms  of  its  deviation  from  c,  and  the  absorp- 
tion coefficient  is  of  order  ± 5 ± 10~14  cm'1  which  is  close  to 
the  analytical  result  0. 


- 
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Table  7.3.  Run  2,  Zero-ir  Pulse  Error  Check 
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0 0 0 0 
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1 . ‘I-',  r .1  ..■.■■13 

6 

6 . 95E+0G 

-6.2F45QSE-0S 
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In  Figures  3a-d  we  show,  as  a function  of  time,  the  pulse 
field,  polarization,  population  and  power  for  the  0 tt  pulse 
considered.  The  frame  is  2 nsec  wide  in  time  and  shifts  from 
one  space  point  to  the  next  to  keep  the  pulse  center  in  the 
center  of  the  frame.  The  points  in  space  cover  up  to  2.67  x 10 
cm  (267  kilometers)  propagation.  Of  note  is  that  the  pulse 
experiences  a bit  more  than  one  period  in  its  cycle  by  the  end 
of  the  distance  propagated.  In  Fig.  3e  we  show  ‘-he  magnitude 
squared  of  the  Fourier  transform  of  the  normalized  electric  field. 
The  frequency  is  measured  in  cm-1.  One  observes  that  the  trans- 
form is  periodic  in  space,  "breathing"  as  it  were  while  propagating. 

In  4a-d  we  show  "three  dimensional"  displays  of  the  nor- 
malized field,  polarization,  population  and  power  for  lossless 
10  psec  Ott  propagation,,  in  time  (T)  and  space  (R)  . 
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Run  3:  Lossy  Otr  pulse  propagation. 

In  this  run  we  use  the  following  parameters  for  the  case 
of  lossy  propagation  of  a Ott  pulse: 

tQ  = 100  psec  X 1 = 6.1  cm”* 

Nq  = 2.77  x 1011  cm"3  Aw  - 0 (7.10) 

U - 4.58  x 10  2 eX  = 185  psec 

t = t^  = 100  psec 

and  for  time  and  space  steps  we  use  the  values  in  (7.2). 

The  primary  results  are  shown  in  Figs.  5a-d.  Note  that 
the  computer  graphics  omitted  the  very  first  values  of  the 
function  plotted.  The  pulse  area  is  initially  close  to  zero 
and  increases  to  slightly  more  than  unity  by  5 x 10®  cm  pro- 
pagation. This  is  due  to  the  absorption  of  the  second  lobe 
of  the  electric  field  which  has  negative  phase.  As  the  pulse 
continues  through  its  cycle  the  area  decreases  and  becomes 

7 

negative  by  1.5  x 10  cm.  The  pulse  energy  is  initially  at 
the  0— tt  value  of  16  in  normalized  units,  and  decreases  to 
about  5%  of  that  value  by  the  end  of  the  run.  In  Fig.  5c 
the  linear  (upper  curve)  and  nonlinear  (lower  curve)  absorption 
coefficients  are  shown.  The  linear  absorption  coefficient  is 
determined  by  solving  first 
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- + — iAui . 

3 


at 


2 j 


PjU,  t)  = - Wj  E (x,  t) 


(7.11) 


for  the  polarization  (this  is  the  linear  limit  of  the  equation 
in  (2.4))  and  then  computing  the  absorption  coefficient  through 
(4.2)  and  dividing  by  (ctQ)  to  get  the  result  in  absolute 
units.  Basically  the  nonlinear  absorption  coefficient  is  ini- 
tially about  0.6  the  value  of  the  linear  equivalent,  and  the 
ratio  increases  towards  unity  as  the  pulse  energy  decreases. 

The  medium  response  becomes  increasingly  linear  as  the  pulse 
is  attenuated.  The  ratio  of  the  absorption  coefficients  is 
shown  in  Fig.  5d. 

In  Figs.  6a-d  we  show  plots  versus  time,  at  different 
spatial  positions,  of  the  field,  population,  power  and  the 
magnitude  squared  of  the  Fourier  transform.  The  electric  field 
is  seen  to  decay,  although  it  manages  to  cycle  through  more  than 
a period  as  in  the  nonlinear  case.  In  6b  the  populations  are 
shown,  demonstrating  that  as  the  pulse  propagates,  the  satura- 
tion decreases  markedly.  The  pulse  power  is  attenuated  as 
shown  in  6c,  and  in  6d  one  can  get  a glimpse  of  what  the  at- 
tenuation means  in  transform  space,  namely  that  the  medium  is 
chewing  out  the  center  of  the  spectrum. 
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Run  4s  Lossy  short  07t-pulse  propagation. 

We  now  consider  propagation  of  a much  shorter  Ott  pulse 
in  the  resonant  absorption  case.  We  assume  the  following  para- 
metric values 

t = 10  psec 

N = 2.77  x 1011  cm' 
o 

V = 4.58  x 10"2  eX 
t = t^  = 10  psec 

and  for  time  and  space  steps  we  use 


min 

0 

10®  cm 

fcmin 

3= 

-130  psec 

max 

3 x 

tmax 

- 

130  psec 

n 

X 

< 

3 x 

10®  cm 

At 

= 

0.52  psec 

nx  = 

100 

steps 

nt 

S 

500  steps. 

In  Figs.  7a-d  the  primary  results  are  summarized.  Fig.  7a 
shows  the  pulse  area  which  makes  excursions  between  0.04  and 
-0.05.  Since  the  pulse  is  shorter  now  than  in  the  previous 
runs,  the  loss  plays  a smaller  role  and  the  pulse  area  stays 
much  closer  to  zero.  The  pulse  energy  is  shown  in  Fig.  7b,  il- 
lustrating near  linear  decay  from  16  to  about  9 in  normalized 

g 

units  during  3 x 10  cm  propagation.  The  linear  (upper  curve) 


X ^ » 6.1  cm  ^ 

Au>  = 0 (7.12) 

= T2  * 185  psec 


and  nonlinear  (lower  curve)  absorption  coefficients  are  shown 
in  7c,  illustrating  a "ringing"  in  the  linear  absorption 
coefficient.  This  result  is  due  to  the  fact  that  the  Fourier 
spectrum  "breathes"  as  it  propagates,  and  the  different  parts 
of  the  cycle  see  different  absorption  coefficients.  The  ratio 
of  the  nonlinear  to  the  linear  absorption  coefficient  is  shown 
in  7d.  In  Fig.  8a-d  we  give  the  field,  population,  power  and 
square  of  the  transform  for  the  run,  as  function  of  time,  at 
different  spatial  positions. 
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Run  5j  Lossy  long  On-pulse  propagation. 

In  this  run  we  example  lossy  long  pulse  propagation,  using 
the  following  parameters 

tQ  = 1 nsec 

N0  = 2.77  x 1011  cm”3 

U = 4.58  x io-2  e& 

te  = tf  = 1 nsec 

and  for  time  and  space  steps  we  use 


X"1  = 6.1  cm'1 

Ao)  - 0 (7.14) 

= Tj  = 185  psec 


x . 
min 

x 

max 

Ax 


0 

3 x io**  cm 

4 

3 x io  cm 


nx  = 100  steps 


t . 
min 

t 

max 


-13  nsec 
13  nsec 


At  * 52  psec 


nt  = 500  steps. 


(7.15) 


The  results  are  shown  in  Fig.  9a-d  in  which  plots  of  the  pulse 
area,  energy,  absorption  coefficients  and  ratio  of  absorption 
coefficients  are  found. 
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VIII.  Results:  One  Line,  Off  Resonant  Case 


Run  1:  Lossless  2 it  pulse  propagation. 


In  the  first  run  we  consider  propagation  of  a 2ir  pulse 
off  resonance,  comparing  analytical  results  with  the  numerical 
results.  In  this  problem  we  take  the  following  parameter  values 


t = 100  psec 
Nq  = 2.77  x 1011  cm"3 
y = 4.58  x 10"^  eR 


X"1  = 6.1  cm"1 
Ait)  - -0.05  cm"1 
T1  = T2  = 1050  sec 


(8.1) 


and  for  time  and  space  steps 


x 

x 


min 

max 
Ax  = 

nx  = 


= 0 


3 x 10  cm 
3 x io^  cm 
100  steps 


(8.2) 


The  input  field  was  2sech  t in  normalized  units.  The  results 
are  given  in  Table  8.1.  The  real  pulse  area  oscillates  sinu- 
soidally, the  pulse  energy  is  within  0.3%  of  the  analytic 
value.  The  pulse  velocity  is  probably  accurate  to  within  1% 
in  the  difference  between  it  and  c.  The  absorption  coefficient 
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Table  8.1 


OLOCH  R * 1.011E+0?  HCYCLE  ■ 301 


VARIABLES  IN  SPACE 


AFEAR  ENERGY  VELOC1TYC 

KAPPA  1 

R 1 

0 0 

0 

0 

2 

2. ICE +06 

5.96188RE+03 

4.63" PIPE 

-87 

6. 289054"- 14 

3 

3.30E+0G 

5.509465E+09 

?.97??7.9£+Q0 

4.%33f:i.  ■ 

-07 

4.067093E - 14 

4 

4.30E+0G 

4.875769E+30 

7 . 9772  i-6E+00 

4. 2333 ! 

-07 

5. 48665!  U 

5 

5 . 762+36 

4.08094431+60 

7 . 977 i i5E+G0 

4 . C i . 

-37 

3.4360 70E-  14 

6 

6.90E+06 

3. 150494E+O0 

7.977 ' K-1+S0 

4. 8330 l; . 

-37 

4.41813-14 

? 

8. 10E+S6 

2.  1 14616E-3  00 

7, 9772  73E+08 

4. 83700r 

-07 

6.  19.365':  - 4 

9 

9.30E+06 

1 .007266E+00 

7.9?i  51^+03 

4.: 

-37 

5. 47200., 7 14 

9 

i . c)rst:-+o?' 

-1.350016E-01 

7.977 : "0E+30 

-87 

4.45450;.  • - i 4 

ID 

1.1 7E+07 

-1.274145E+00 

7.977?  r„)E+33 

4.830001: 

-07 

5. 3 91277c  -’  4 

1 l 

1.29E+07 

-2.372004E+00 

7.9772  18E+0& 

4.832991 

-07 

C n s*"  • : -t  : > 

I TJ 

1.41E+07  . 

-3.3.91639E+00 

7.377-.  J7E+30 

4.93299?- 

-97? 

:>•*,!  C.  '.4  . . -1 

! .53C+07 

-4 . 298G72E+00 

7.3772  17E+30 

4.  11 

-07 

- / * 4 . lit"  : ••( 

14 

1 . 651-  +0? 

-5. 062542E +00 

7.377  108+33 

4,.  @3303 1 

-97 

3.826252E-14 

1 13 

1 . ?r’E+0r’ 

-5.657S96E+00 

7.977  2 :;5E+00 

4 . 53"  7 

-0? 

4.074611;  ’< 

16 

1 .R9I.+0? 

-6.063974E+00 

7.3 77  - -3E+00 

4.232::?:. 

-07 

4.c  9<  ■ 

1 7 

2.O1E+07 

-6.268262E+80 

7.377  J SE+00 

4.32-001 

-07 

3.41  6221 E- . ; 

la 

2. 13E+07 

-6 . 26381S3E+S3 

7.97.72  I5E+00 

4.-3303?' 

1 . . r63 . . 

* IT) 

2.25E+3? 

-6.051731E+B0 

7 . 97  7 K3E+O0 

4 933007: 

-37 

7.?' - 0575  >14 

2U 

2.U7E+Q7 

-5 . 637773E+00 

7.9772  iSE+00 

4.533S'?' 

-37 

2.81743'  -.4 

21 

2.49E+8? 

-5 . 037S67E+00 

7.977  ! 15E-  23 

4.33301 

37 

2.0  7: 99?  4 

72 

•?.fi  1E+07 

-4.2  7Q845E +30 

7.977  3 OE+O0 

4.43re-;r‘ 

-07 

2.62477  - 

A3 

2.73E+0? 

-T.3fi3??tE+©3 

6E  1-09 

4.8333  ’ ' 

J? 

2.49G8 ; : 

24 

.831  t07 

-2.3430U5E+00 

;.p#0 

4.  ■ 

...  j? 

2.02571- 

7S 

2. y vt +0? 

-1.250029E+00 

7.377,.  .fiErQb 

4.83 

• J? 

2.8721* 

/ 

/ 


is  less  than  10 
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cm  The  accuracy  of  the  calculation  is 
sufficient  for  the  purposes  of  the  analysis. 

Off-Resonant  Propagation  of  "Zero-pi  Pulse" 

As  explained  in  the  appendix,  there  is  no  zero-pi  pulse 
like  solution  for  propagation  in  a collision  free  medium,  i.e.  a 
pulse  with  its  spectrum  not  centered  at  the  material  line. 
Computer  runs  made  with  a zero-pi  pulse  with  a spectrum  com- 

I 

puted  at  resonance  and  then  simply  shifted  off  resonance  show 
break-up  of  the  pulse  into  two  separate  pulses. 

When  loss  is  introduced,  the  pulse  does  not  necessarily 
break  before  linear  propagation  is  reached  due  to  the  loss. 

In  this  set  of  runs  we  consider  the  case  of  lossy  propagation 
at  several  values  of  the  detuning  Aw.  In  lOa-d  we  show  the 
Fourier  transform  magnitude  squared  of  the  normalized  field  in 
space  (cm)  and  frequency  (cm-1),  at  Aw  * 0,  .05  cm-1, 

.10  cm"^  and  .20  cm  Basically,  when  the  pulse  is  chirped, 
one  of  the  lobes  (this  one  which  has  more  overlap  with  the  lossy 
line)  is  eaten  away  faster  than  the  other. 

In  Fig.  lla-d  are  shown  the  normalized  field,  polarization, 
population  and  absolute  power  for  a chirped  Ott  pulse  defined 
to  Aw  = 0.05  cm  1.  Of  note  is  that  the  pulse  does  not  break 
up — a result  which  follows  simply  from  the  fact  that  breaking 
up  is  a nonlinear  effect,  and  the  pulse  does  not  remain  nonlinear 
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long  enough  to  break  up  (te  = = 100  psec) . In  Fig.  12a-d 

are  presented  the  normalized  electric  field  versus  T and  R 
at  Aw  = 0,  .10  cm  \ .20  cm  \ 0.60  cm""^ . One  can  see  that 
the  nonlinear  gyrations  of  the  field  are  slower  in  the  case  of 
detuned  propagation;  when  the  line  is  detuned,  the  effect  of 
the  line  on  the  pulse  becomes  less. 

In  13a-c  are  displayed  the  normalized  real  polarizations 
(Imaginary  polarizations  are  not  shown).  The  effect  of  the 
field  on  the  medium  is  less  off  resonance.  In  Fig.  14  are  shown 
the  populations  (normalized)  at  Au)  = 0,  0.1  and  0.2  cm-1. 

The  saturation  is  much  less  off  resonance,  and  the  "rippling" 
which  is  a nonlinear  effect  is  virtually  absent  off  resonance. 

Fig.  15  shows  the  pulse  power  at  Aw  = 0,  .10  cm-i  and 

.20  cm  One  observes  simply  that  the  rippling  is  less  at 
resonance.  The  message  is  that  as  one  gets  further  off- 
resonance,  the  0— ir-like  nonlinear  effect  (i.e.  rippling)  is 
weaker.  The  difference  between  the  linear  absorption  and  non- 
linear absorption  is  primarily  a saturation  effect  rather  than 
due  to  any  other  special  Oir  property. 

Figs.  16,  17  and  18  give  plots  versus  distance  respectively 
at  Aw  = 0,  0.05,  and  0.10,  for  tg  = tf  = 100  psec. 

(a)  pulse  area  (real  and  imaginary) 

(b)  normalized  pulse  energy 

(c)  the  deviation  of  the  pulse  velocity  u from  the 
speed  of  light  in  terms  of  1 - u/c 
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<d) 


(e) 


the  loss  coefficients  for  linear  and  nonlinear  absorption 
the  ratio  of  nonlinear  to  linear  loss. 


IX.  H^Q  Lines  Adjacent  to  CH^F  and  D^O  Laser  Lines 


The  Methyl  Fluoride  (CH^F)  laser  at  496  u and  the 
D20  laser  at  385  p are  two  available  submillimeter  sources 
whose  radiations  fall  between  main  absorption  lines  of  H20. 
For  this  reason  we  study  specifically  the  two  radiations  and 
represent  the  water  vapor  absorption  by  the  nearest  adjacent 
lines. 

Water  is  an  "asymmetric-top"  molecule  (Townes-Schawlow, 
1955) . The  energy  levels  of  such  a molecule  are  described  in 
terms  of  three  numbers,  J,  K_^  and  respectively.  J 

is  the  total  angular  momentum  and  is  a "good"  quantum  number. 
K_^  and  are  not  "good"  quantum  numbers  in  that  they  do 

not  represent  the  angular  momentum  around  one  of  the  three 
principal  axes  of  the  molecule.  In  general,  there  is  no  level 
degeneracy  for  a particular  set  of  values  of  J,  K_^»  K!' 

There  is,  however,  an  orientational  degeneracy. 

We  obtained  line  strengths  from  the  computer  printouts  of 
S.  A.  Clough.  These  line  strengths  contain  information  on  the 
dipole  moment  <|p|2>  averaged  over  all  orientations  of  the 
molecules.  One  may  start  with  the  two  level  system  equations, 
(1.15)  and  (1.16)  of  the  Appendix,  in  the  small  signal  limit 
in  which  - p22  can  be  treated  as  constant.  Then,  one 
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solves  for  p12  on  line  center  and  puts  the  result  into  (1.10) 
of  the  Appendix.  The  spatial  rate  of  change  of  E may  be 
equated  to  the  amplitude  decay  coefficient  (a  in  the  notation 
of  Marcuse,  1)70,  p.  286,  not  to  be  confused  with  the  a used 
here  earlier).  One  finds,  in  agreement  with  Marcuse, 


a 


uj  | u | 2 T2 
2 ech 


N(p 


22 


pn) 


(9.1) 


An  orientational  degeneracy  forces  one  to  take  averages  over 
all  M so  that  the  average  a,  <a>,  contains  the  average  of 

M2,  <M2>. 

The  line  strengths  s,  as  used  by  Clough,  are  the  <a> 
integrated  over  the  line  width  function 


g (v) 


1 + 


1 


Av 


(9.2) 


where  2ttAv  = 1/T2»  and  expressed  in  inverse  centimeters. 
The  conversion  factor  is 


1 sec"1  = 3.34 


x lO-’11  cm-1. 


(9.3) 


Combining  all  this  one  finds 


■■■■■■ 
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uD  = /<|u|2>  = 


sX 

N , 
abs 

N„  « 

7.85  x 10-15 

1/2 


debye 


(9.4) 


with  NabS/NH20  “ p22  ~ P11  defining  the  ratio  of  the  popula- 
tion density  difference  participating  in  the  transition  to  the 
total  population  density  of  water  molecules. 

Further, 


WV  - «p  - 


hv. 


kT 


1 - exp  - 


hv  -hv. 
u l 


kT 


(9.5) 


Any  molecule  of  angular  momentum  J can  have  2J  + 1 orien- 
tations with  respect  to  a fixed  spatial  axis  (say  the  applied 
electric  field  direction) . These  orientations  are  described 
by  the  "good"  quantum  number  M where  Z is  the  partition 
function.  For  H20  at  T = 296°K,  Z = 174.83  (S.  A.  Clough, 

private  communication) . 

The  constant  g is  the  statistical  weight  of  the  lower 
level  ( Townes -Schawlow,  Table  4-7) 


g = 1 for  symmetric  level 
g = 3 for  antisymmetric  level. 
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The  symmetry  or  antisymmetry  is  determined  from  K_ ^ and  Kx , 
depending  upon  whether  K-;l  + K1  is  even  or  odd.  Further, 
hvQ  = hvu  - hv^,  where  hvu  is  the  energy  of  the  upper  level, 
hv^  that  of  the  lower  level.  Now,  <|u|2>  gives  the  average 
matrix  element  for  absorption  of  a field,  with  fixed  orientation 
with  respect  to  a spatial  axis,  say  the  z-axis.  The  z-component 
of  u has  two  different  dependences  upon  M,  one  for  a transi- 
tion, J J'  <=  J + 1,  the  other  for  J -*■  J'  = J.  in  the 
Fourier  case  (Townes-Schawlow,  1-73) 

|u  |2  - (J  + l)2  - M2;  J'  = J + 1 (9.6) 

Z 

In  the  latter  case, 

|yzl2  « M2;  J’  = J.  (9.7) 


We  find  by  simple  recursion 


E M2  = — J (2J  + 1) (J  + 1) . (9.8) 

-J  3 


The  initial  population  difference  for  nonlinear  interaction 
with  any  one  of  the  transitions  of  given  M is 


N 


h2o 


g(2J  + 1) 


exp  - 


h V( 


kT 


1 - exp  - 


fhvu  - hVjl) 


kT 


1* 

g.v. 


where 


gA/gu  = 1 for  AJ  = 0 

9t/9„  = — -*-1  for  AJ  =*  1. 

* U 2J  + 3 

The  reason  for  this  is  that,  in  the  absence  of  collisions, 
a transition  J -*■  J + 1,  M M,  reaches  only  the  M'  = M 
states,  and  not  the  final  ones  M’  = + (2J  + 3).  These  states 
are  unaffected  by  the  radiation  and  get  filled  only  after  M- 
changing  collisions  have  occurred.  Such  collisions  cannot  be 
dipole-dipole  collisions,  and  thus  they  occur  at  a much  lower 
rate  than  1/T2.  We  shall  ignore  these  collisions  and  hence 
ignore  the  M'  = ± (2J  + 3)  states  of  the  upper  level  in  our 
propagation  studies. 

The  parameters  used  in  the  computations  are  summarized  in 
Tables  9.1  and  9.2.  We  studied  resonant  excitation  of  the  four 
lines  at  18.577,  20.704,  25.080  and  30.560  cm  \ that 

lie  near  the  CH3F  and  D20  laser  lines.  The  new  feature 
is  that  these  lines  are  degenerate  and  therefore  the  different 
M-levels  are  excited  differently.  There  is  no  zero-pi  pulse 
solution  in  the  loss-free  case,  strictly  speaking.  We  picked 
a pulse  that  would  be  a zero-pi  pulse  in  a nondegenerate  system 
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Table  9.1 


Summary  of 

H2° 

Data  on 

4 Lines 

Line 

<JK-lKl>lo»er 

(JK  -K.) 

-1  1 upper 

g 

F* 

F 

u 

nl 

E 

..  “ kT. 

U e l 

Nabs 

-1 

cm 

-1 

cm 

-1 

cm 

x 1015 

-2 

x 10 

x 101* 

18.577 

110 

101 

3 

23.794 

42.371 

9.04 

8.62 

7.79 

20.704 

532 

441 

3 

488.1 

508.1 

3.48 

9.56 

3.33 

25.080 

211 

202 

1 

70.091 

95.171 

4.01 

11.5 

4.61 

30.560 

422 

331 

1 

285.219 

315.779 

2.54 

13.8 

3.51 

Population  density  of  lower  level: 

-E/kT 


N. 


gN(2J  + 1)  e 


Z - 174.83 


Population  difference: 


abs 


Nl(1  - e 


-AE/kT 


) 


Line 

N . /Nu  . 
abs  1^0 

M 

s(n  . : 

abs 

-1 

cm 

per  molecule  x cm 

-3 

x 10 

-2 

debye  x 10 

per  molecule  x ( 

18.577 

-20 

5.00  x 10 

3.95 

2.95 

1.27 

20.704 

-22 

5.65  x 10 

1.69 

0.454 

0.0334 

25.080 

-20 

3.47  x 10 

2.34 

2.74 

1.48 

30.560 

-21 

1.43  x 10 

1.78 

0.579 

0.0803 

The  relaxation  time: 


T2(sec) 


3.34  x IQ 
2Av 


-11 


Av  in  cm 


-1 
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Table  9.2 

Summary  of  Input  Data  for  Individual  Levels 

Case  1;  J = 1,  J'  = 1. 


Line 

M 

Nabs  x densifcy 

T2 

-1 

cm 

debye  x 10  1 

_ -3  ,.14 

cm  ^ 10 

sec  x io  '*'2 

18.577 

2.95 

7.79 

155 

1 M | = 1 

4.43 

5.19 

Case  2 : J = 

5,  J'  = 4 

Line 

Nabs  x densifcy 

T2 

cm  ^ 

debye  x 10~2 

cm-3  x io13 

-12 

sec  x 10 

20.704 

4.54 

33.3 

230 

1 M | = 3 

3.96 

7.40 

1 M | = 2 

5.20 

7.40 

3 

II 

H* 

5.94 

7.40 

| 


I M | = 0 


6.19 


3.70 


with  the  same  value  of  /<|y|2>.  Note  that  the  power  in  such 
a pulse  scales  like  l/< | u | 2>tez • Because  the  <|y|2>  values 
of  the  lines  are  not  greatly  different,  we  are  dealing  in  all 
cases  with  pulses  of  peak  power  of  the  order  of  10  Mw/cm2  for 
100  psec  pulse-length — comparable  to  the  pulse  power  of  the 
6.1  cm  1 line  analyzed  earlier. 

The  18.577  line  does  not  represent  much  novelty.  It 
behaves  like  a nondegenerate  line  because  only  the  levels 
M =1  interact  with  the  radiation  and  they  have  the  same 
dipole  moment.  The  only  feature  worth  noting  is  that  the  pene- 
tration dipole  is  very  short.  The  distance  scale  terminates  at 
10  m — the  absorption  line  is  indeed  very  strong.  Saturation 
has  reduced  the  loss  to  55%  of  its  nonlinear  value,  but  after 
about  6 m propagation  the  loss  has  approached  the  linear  loss 
to  within  85% — the  penetration  has  been  improved,  but  very 
little. 

Figure  21  plotted  for  the  20.704  cm  ^ line  is  on  a scale 
of  1.4  km — a much  weaker  absorption.  We  note  that  the  dif- 
ferent M-levels  are  acted  upon  differently  by  the  radiation, 
and  the  loss  in  each  is  reduced  by  different  amounts  due  to 
nonlinear  effects.  The  transition  is  one  with  AJ  = 1 so  that 
the  strongest  dipole  moment  is  for  M = 0.  Figure  21d  shows 
that  indeed  the  largest  reduction  in  loss  occurs  for  M = 0, 
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to  about  45?.  of  the  linear  value.  For  the  M = 0 line  the 
equivalent  "average"  zero  pi  pulse  used  for  the  computation  has 
more  energy  than  a true  zero-n-pulse  would  have  for  this  level 
of  strongest  dipole  moment,  hence  the  large  reduction  in  loss. 

Yet,  after  roughly  a 1/e  decay  length  the  loss  reduction  has 
disappeared  for  all  practical  purposes. 

Figure  22  is  for  a line  with  AJ  = 0.  The  strongest  dipole 
moment  is  for  |m(  = 2.  From  Fig.  22c  we  gather  that  indeed 
this  level  has  the  larger  loss  reduction,  to  about  60%  of  its 
linear  value.  Note  that  the  penetration  depth  is  of  the  order 
of  5 m,  again  this  is  a very  strong  absorption  line. 

Figure  23  is  for  a AJ  = 1 transition  and  qualitatively 
duplicates  the  behavior  of  the  AJ  = 1 transition  at  20.704  cm-*. 
Figure  21.  It  is  a weaker  absorption  line,  the  horizontal  scale 


extends  to 


450  m. 
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X.  Propagation  of  CH^F  Laser  Radiation 

In  the  preceding  section  we  studied  propagation  of  radiation 
in  resonance  with  the  degenerate  lines  adjacent  to  the  CH.jF  and 
D2O  lines  at  495  y and  385  y respectively.  The  propagation 
characteristics  were  not  greatly  different  from  those  in  resonance 
with  a nondegenerate  line.  As  a final  test,  we  study  the  propaga- 
tion of  495  y radiation  as  it  is  affected  by  absorption  via  the 
strongest  adjacent  line  , the  18.577  cm-1  line.  Figs.  24a, b 
show  the  propagation  of  a 10  psec  pulse.  The  shape  of  the 
pulse  is  taken  as  that  of  a zero-pi  pulse  on  resonance  with  the 
line,  detuned  to  495  y = 20.2  cm”*.  Curve  A of  Fig.  24a  has 
the  energy  appropriate  for  a zero  pi  pulse,  the  others  are  of 
higher  energy.  The  1/e  decay  is  roughly  500  m.  Fig.  24b 
shows  several  interesting  features.  Curve  A is  for  a small 
signal  pulse.  It  shows  a dramatic  decrease  of  k — the  attenua- 
tion constant,  as  a function  of  propagation.  This  is  due  to  the 
fact  that  the  portion  of  the  spectrum  that  lies  on  the  high 
absorption  side  of  the  line  gets  "eaten-up"  initially.  The 
remaining  part  of  the  spectrum  has  smaller  loss  as  it  continues 
its  propagation. 

I 

Curve  B of  Fig.  24b  corresponds  to  curve  A of  Fig.  24a. 

The  initial  saturation  effect  is  very  noticeable  reducing  k to 
about  30%  of  its  small  signal  value.  As  the  pulse  propagates, 

I 

J 
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the  absorption  coefficient  exceeds  the  small  signal  value 
because  the  most  "favored"  part  of  the  spectrum,  nearest  to  the 
line  center  of  the  absorption  line,  has  a weaker  saturation  ef- 
fect after  attenuation  has  occurred. 


Conclusions 


The  zero-pi  pulse  as  published  by  Lamb  (1974)  is  a special 
case  of  a self  induced  transparency  solution  obtained  when  the 
spectral  line  of  the  medium  is  symmetric  and  the  carrier  fre- 
quency of  the  pulse  is  centered  with  respect  to  the  line.  Fol- 
lowing "small  signal"  reasoning,  we  were  interested  in  detuning 
pulses  from  line  center  so  that  the  loss  is  avoided  at  least 
partially  when  the  "center  of  gravity"  of  the  spectrum  lies 
outside  the  absorption  line.  Generalized  solutions  have  not 
been  published  although  the  formalism  of  inverse  scattering 
theory  can  be  applied  to  that  case.  Lamb,  in  another  paper 
(1973),  gave  a general  formula  which  is  wrong  (his  equation  16). 
For  this  reason  the  analysis  was  redone  to  provide  the  starting 
point  for  the  computations  presented  in  this  report. 

We  have,  therefore,  applied  the  Inverse  Scattering  Theory 
to  obtain  an  expression  for  the  generalized  zero-pi  pulse.  In 
the  course  of  this  work  we  found  it  necessary  to  cast  the 
formalism  developed  by  mathematicians  into  language  and  models 
understandable  to  the  physicist  and  engineer.  A paper  on  the 
first  version  of  this  work  will  appear  in  Reviews  of  Modern 
Physics.  The  appendix  of  this  report  is  a further  simplifica- 
tion which  avoids  entirely  the  use  of  the  Marchenko  equation 
of  Fourier  transform  space  and  thus  gives  more  rapidly  results 
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useful  to  the  evaluation  of  the  zero-pi  pulse.  One  important 
outcome  of  the  analysis  is  that  there  exists  no  zero-pi  pulse 
detuned  from  line  center.  Any  pulse  of  the  "breather  type" 

(such  as  the  zero-pi  pulse  is  called  because  it  changes  shape 
as  it  propagates)  breaks  up  when  its  spectrum  is  not  centered 
with  the  line  of  the  medium.  Therefore , we  used  in  the  computer 
analysis  in  the  off-resonant  case  an  electric  field  envelope 
evaluated  from  the  on-resonant  case,  and  shifted  off  resonance. 
Such  pulses  are  found  not  to  break  up  when  sufficient  loss  is 
introduced  so  that  the  intensity  of  the  pulse  is  decreased  to 
a sufficiently  low  level  before  the  break-up  had  a chance  to 
occur . 

The  numerical  studies  proceeded  in  two  stages.  First,  we 
looked  at  2tt  and  O-ir  pulse  propagation  in  the  absence  and 
presence  of  loss  at  and  near  the  6.1  cm-1  line  of  H20,  with- 
out including  the  degeneracy  of  the  level  so  that  ideal  2tt 
pulse  and  0-tt  pulse  propagation  in  the  absence  of  loss  was 
possible.  This  numerical  study  showed  certain  features  of 
zero-pi  pulse  propagation  in  the  presence  of  loss  in  the  simplest 
possible  context.  The  analysis  also  enabled  one  to  test  the 
program  for  its  potential  of  predicting  long  distance  propaga- 
tion (up  to  3 x 107  cm) , in  this  case  of  low  absorption.  The 
following  specific  findings  are  of  interest,  but  note  that 
level  degeneracy  had  not  been  taken  into  account — thus  the 
beneficial  effects  of  saturation  are  exaggerated. 
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(a)  The  attenuation  of  the  6.1  cm  ^ 1^0  line  is  reduced 

to  .6  of  its  linear  value  for  a zero-pi  pulse  over 

a distance  within  which  the  energy  has  not  been  greatly 
depleted.  After  that  the  pulse  propagates  in  essentially 
a linear  manner.  For  a 100  psec  pulse,  the  required  peak 
intensity  is  8 MW/cm2.  Typical  propagation  distances  are 
20  km. 

(b)  The  off-resonant  "zero-pi  pulse"  can  lead  to  a further 
reduction  of  loss  to  .25  of  its  linear  value,  if  one 
peak  of  the  double-peaked  spectrum  of  the  zero-pi  pulse  is 
made  to  coincide  with  the  absorption  line.  This  is  an  in- 
dication of  a saturation-bleaching  effect,  the  bleaching 
is  more  effective  when  the  spectrum  has  greater  net 
intensity  within  the  spectral  response  width  of  the 
medium. 

(c)  Further  detuning  of  the  "zero-pi"  pulse  spectrum  from 
the  absorption  line  gives  less  loss  reduction. 

(d)  As  a pulse  is  detuned  from  line  center,  the  loss  is 
reduced  as  one  would  expect  from  a linear  analysis  be- 
cause less  of  the  spectrum  coincides  with  the  absorption 
line.  In  fact,  a "zero  pi"  pulse  of  8 MW/cm2  intensity 
as  defined  for  on  resonance  propagation,  when  detuned  by 
more  than  a linewidth,  interacts  with  the  medium  in  an 
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almost  linear  fashion.  This  is  because  a relatively  small 
portion  of  the  spectrum  interacts  with  the  medium. 

(e)  The  zero-pi  pulse  has  no  small  signal  limit.  Accordingly, 
diffraction  effects  will  have  an  effect  on  zero-pi  pulse 
propagation  within  a beam  of  finite  cross  section  similar 
to  that  of  2 tt  pulses. 

(f)  The  CH^F  and  laser  lines  lie  in  the  wings  of  very 

strong  water  vapor  absorption  lines.  These  lines  connect 
degenerate  levels  and  hence  S.I.T.  solutions  do  not  exist 
in  the  lossless  (collision-free)  case  (except  for  the 
18.577  line).  Using  on  resonant  propagation  at  any  one 
of  the  four  lines  of  "equivalent"  zero-pi  pulses  with 

(+  -)  areas  chosen  to  correspond  to  the  average  matrix 
element  of  the  transition,  one  observes  reductions  of 
absorption  due  to  saturation  very  similar  to  those  of  the 
nondegenerate  case,  of  the  order  of  55%  of  the  linear 
loss.  The  penetration  depths  are  now  much  shorter;  of 
the  order  of  5m  for  the  two  strong  lines,  500  m for 

v 

the  two  weak  lines. 

(g)  The  peak  power  of  the  equivalent  zero-it  pulse  scales  like 
1/ | m | 2 tQ2.  Hence  the  peak  power  one  needs  to  experience 
nonlinear  loss  reduction  is  of  the  order  of  10  Mw/cm2  for 
all  the  lines — their  dipole  moments  being  all  of  order 

0 . 5 debye . 
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(h)  The  and  D2O  laser  lines  lie  in  a regime  of 

relatively  weak  absorption.  Fig.  1 shows  that  the  linear 
loss  exponent  2alQ  within  1 km  is  9.5  and  16  re- 
spectively, corresponding  to  1/e  decay  distances  of 
105  m and  62  m respectively. 

(i)  On-resonant  propagation  on  the  20.704  cm*1  line  produces 
1/e  decay  distances  of  the  order  of  360  m which  is 
much  greater  than  the  105  m distance  of  linear  propaga- 
tion read  off  the  graph.  The  discrepancy  is  due  to  the 
fact  that  the  wing  of  the  strong  adjacent  line  at  18.577 
influences  the  linear  loss  at  the  20.02  cm  1 line  of 
the  CH^F  laser  reducing  the  1/e  decay  distance  ap- 
preciably. 

(j)  Off-resonant  propagation,  at  the  frequency  of  the  CH^F 
line,  of  a zero-pi-like  pulse  of  10  psec  duration  shows 
the  following  characteristics:  In  the  small  signal  limit, 
the  spectral  component  of  the  pulse  on  the  high  side  of 
the  absorption  line  gets  eaten  up  and  the  attenuation 
coefficient  decreases  with  further  propagation.  An  in- 
tense pulse  achieves  appreciable  initial  reduction  of 
attenuation  (to  30%  of  small  signal  limit)  until  the 
absorption  has  reduced  its  intensity. 
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(k)  In  the  case  when  the  pulse  energy  is  several  times  the 
on-resonant  zero-pi-pulse  energy,  the  pulse  breaks  up 
into  a series  of  pulses  traveling  at  different  speeds. 

(l)  The  zero-degree  pulse  propagation  with  low  loss  as  observed 
by  Grieneisen  et  al.  was  a linear  phenomenon,  associated 
with  the  peculiar  spectral  distribution  characteristic  of 

a zero  degree,  small  signal  pulse  (which  is  also  character- 
istic of  the  zero-pi  large  signal  pulse) . The  spectrum 
has  zero  intensity  at  the  carrier  frequency.  Hence  a zero- 
degree  pulse  with  the  carrier  centered  at  the  line  exper- 
iences reduced  absorption  because  the  "bulk"  of  the  spectrum 
lies  outside  the  absorption  line.  There  are  other  subtle 
effects  in  the  limit  of  inhomogeneous  broadening,  that 
were  discussed  by  Grieneisen  et  al.,  but  do  not  concern 
us  here  because  H2O  is  homogeneously  broadened  at  atmo- 
spheric pressure. 

(m)  The  overall  conclusion  is  that  lower  loss  can  be  achieved 
in  the  linear  regime  when  the  spectrum  is  moved  away  from 
the  frequencies  of  maximum  absorption.  The  use  of  carrier 
frequencies  that  are  far  from  the  peak  of  the  absorption 
lines  are  equally  or  more  effective.  The  zero-degree 
pulse  would  be  useful  only  in  those  cases  in  which  the 
use  of  a laser  frequency  coincident  with  an  atmospheric 
absorption  line  were  required,  e.g.  when  a CC>2  laser 
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is  used  and  absorption  by  atmospheric  C02  is  to  be 
reduced . 
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Appendix  I 

Derivation  of  Two-pi -pulse  and  Zero-pi-pulse 
by  the  Inverse  Scattering  Method 


Abstract; 

The  Inverse  Scattering  Method  (ISM)  as  applied  to  Self- 
Induced  Transparency  (SIT)  is  reviewed.  It  is  shown  that  the 
(Zakharov-Shabat)  scattering  equations  follow  directly  from 
the  two  level  system  equations.  The  analytic  continuation  of 
the  Zakharov  Shabat  equations  into  the  complex  plane  is  inter- 
preted by  an  analog;  the  excitation  of  a parametric  oscillator 
by  sources  growing  exponentially  in  time.  This  model  gives  a 
"physical  feel"  for  the  behavior  of  the  eigenfunctions  of  the 
scattering  problem.  The  spatial  invariance  of  the  eigenvalues 
is  proven  following  a method  outlined  by  Ablamovitz  et  al.  (1974). 
The  2tt  pulse  and  zero-pi  pulse  follow  immediately  and  concisely 
from  the  preceding  analysis.  It  is  shown  that  no  zero-pi -like 
solution  exists  with  a spectrum  centered  off  the  material  line. 
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Introduction 

The  inverse  scattering  method  (I.S.M.)  is  one  of  the  very  few 
general  methods  of  solution  of  a class  of  nonlinear  differential 
equations  in  one  spatial  dimension  and  time  (Gardner  et  al.,  1967; 
Zakharov  and  Shabat , 1972;  Whitham,  1974;  Ablowitz  et  al.,  1974a; 

Lamb,  1973).  In  the  application  of  the  I.S.M.  a linear  scattering 
problem  of  quantum-mechanical  nature  is  associated  with  the  nonlinear 
differential  equation  (Whitnam,  1974).  The  sought-for  solution  of 
the  nonlinear  differential  equation  at  the  initial  time  t = 0 plays 
the  role  of  the  scattering  potential  or  "well”  of  the  linear  scat- 
tering problem.  The  initial  conditions  of  the  problem  to  be  solved 
prescribe  the  transmission  and  reflection  coefficients  of  the  linear 
scattering  problem  from  which  the  scattering  "well"  may  be  determined 
by  standard  techniques  of  (inverse)  scattering  theory.  The  evolution 
in  time  of  the  scattering  problem,  which  may  take  several  forms 
(Ablowitz  et  al. , 1974a, b),  then  prescribes  the  evolution  in  time 
of  the  solution  of  the  nonlinear  equation — or  alternately — the  form 
of  the  nonlinear  differential  equation  associated  with  this  particular 
scattering  problem  (see  Fig.  2). 

In  1973  Lamb  showed  (Lamb,  1973)  how  one  may  associate  with  the 
equations  of  self-induced  transparency  (S.I.T.)  one  of  the  standard 
equations  of  the  I.S.M.,  the  Zakharov-Shabat  equations  (Zakharov  and 
Shabat,  1972).  He  proceeded  through  a set  of  variable  transformations 
with  no  apparent  physical  interpretation.  In  fact,  one  of  the  in- 
triguing unsolved  problems  of  the  I.S.M.  is  the  development  of  a 


70 


procedure  to  find  the  scattering  problem  associated  with  a particular 
differential  equation. 

In  many  cases  the  scattering  problem  of  quantum-mechanical  char- 
acter will  have  no  physical  interpretation.  It  is  to  be  expected, 
however,  that  in  those  cases  in  which  the  physics  underlying  the 
nonlinear  differential  equation  is  based  on  quantum  mechanics  the 
associated  scattering  problem  must  have  a direct  physical  meaning. 

In  the  case  of  the  Josephson  transmission  line  and  S.I.T.,  this  has 
been  pointed  out  by  McLaughlin  and  Corones  (McLaughlin  and  Corones , 
1974).  The  physical  interpretation  of  the  scattering  problem  ob- 
tained after  Lamb  had  reduced  the  S.I.T.  problem  to  one  amenable  to 
the  inverse  scattering  method  is  but  one  example  of  various  develop- 
ments presented  by  different  authors  in  the  course  of  time.  The 
present  author  attempted  to  gain  an  understanding  of  the  I.S.M.  by 
uncovering  physical  interpretations  for  the  mathematical  steps.  Even 
though  most  of  the  specifics  have  appeared  in  the  literature  (Gardner 
et  al.,  1967;  Zakharov  and  Shabat,  1972;  Whitham,  1974;  Ablowitz  et 
al.,  1974a, b;  Lamb,  1973;  McLaughlin  and  Corones,  1974;  Ablowitz  et 
al.  1973;  Faddeyev,  1962),  this  paper  may  serve  as  a guide  for  the 
physically  inclined  to  an  important  mathematical  method. 

The  Zakharov-Shabat  equations  written  as  differential  equations 
in  normalized  time  x with  the  electric  field  £(x,  z)  as  the 
scattering  "well"  define  the  linear  scattering  problem  associated 
with  S.I.T.  The  spatial  variable  z plays  the  role  of  a parameter. 
In  Section  I we  show  the  direct  connection  between  the  Zakharov- 
Shabat.  equations  and  the  equations  of  the  two-level  system  ini  or- 
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acting  with  the  electric  field  in  S.I.T. 

The  nature  of  the  solutions  to  the  Zakhar ov-Shabat  equations 
may  be  anticipated  by  recognizing  their  similarity  with  a well 
known  physical  problem:  the  parametric  interaction  of  two  waves 
via  a nonlinear  medium  excited  by  a pump  wave  (Yariv,  1976).  This 
is  done  in  Section  II.  In  particular,  it  is  known  that  "unstable" 
solutions  growing  in  time  are  encountered  in  parametric  interactions. 
This  fact  may  be  used  to  predict  the  location  of  the  eigenvalues 
of  the  Zakharov  Shaba t equations  in  the  complex  plane. 

The  exploration  in  Section  III  of  the  behavior  of  the  Zakharov- 
Shabat  equation  in  the  complex  plane  leads  directly  to  an  equation 
for  those  solutions  of  the  Zakharov-Shabat  equations  that  correspond 
to  a reflection-free  well,  and  the  shape  of  the  well. 

In  Section  IV  we  show  that  the  assumption  of  independence  of 
z of  the  eigenvalues  of  the  Zakharov-Shabat  equations  leads  to  a 
form  of  the  Maxwell-Bloch  equations  of  S.I.T.  In  Section  V we 
27T-pulse  ar.d  zero-ir  pulse  solutions. 


obtain  the 
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I.  The  Zakharov-Shabat  Equations  as  the  Equations  of  the  Two- 
Level  Systems 

The  Zakharov-Shabat  equations  are  central  to  the  inverse 
scattering  method  applied  to  self  induced  transparency.  In  this 
section  we  review  briefly  the  equations  of  a two  level  system 
excited  by  an  E-field  and  show  that  the  resulting  equations  are 
equivalent  to  the  Zakharov-Shabat  equations  (Zakharov  and  Shabat, 
1972)  arrived  at  by  Lamb  (Lamb,  1973). 

In  the  slow  envelope  approximation,  the  wave  equation  for  the 
electric  field  envelope  E(x,  t)  of  a plane  wave  propagating  in 
the  x-direction  is  in  mks  units  (compare  (Lamb,  1973)) 

3l  + k5l.i“t  (i.i) 

3x  c 9t  2 c e 

where  u>  is  the  "carrier"  frequency,  c the  speed  of  light,  c 
the  dielectric  constant,  and  P the  polarization  in  the  medium. 

P and  E are  parallel  to  each  other  and  transverse  to  the  x- 
direction.  The  polarization  of  the  medium  is  obtained  from  the 
analysis  of  two  level  systems  with  a distribution  of  energy-level 
spacings.  Denote  the  amplitude  of  the  wave  function  of  the  upper 
level  (1)  by  a^,  that  of  the  lower  level  (2)  by  a 2>  One  may 
write  down  two  differential  equations  for  the  amplitudes  and 

a2  as  coupled  by  the  E-field  (Vuylstekc,  1960).  Factoring  out 
the  natural  time  dependences  and  retaining  only  the  slowly  time 

varying  portions  of  the  variables,  one  huF 
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• i • - x5t  /-l  <\ \ 

al  = - B • Pl2  e a2  (1.2) 

ft 

*2  - z E*  • 5!2  e'iJt  *1  (1-3) 

ft 

where  6 = u>  - <d^2,  and  p^2  t^le  matrix  element  between  the 

two  levels.  If  we  define 


(1.4) 


(1.5) 


we  obtain 


dv 


1 


+ 


dt 


E 


12 


ft 


(1.6) 


(1.7) 


These  are  already  the  Zakharov-Shabat  equations,  except  for  a 
normalization.  The  density  matrix  p is  the  statistical  average 
of  the  products  of  the  amplitudes  a^ , a2  or  v^ , v2  and  their 


complex  conjugates.  Since  we  are  dealing  here  with  a pure  state 
(no  collisions)  no  statistical  average  need  be  performed. 


pij  = vi  V*  u-8) 

The  positive  frequency  portion  of  the  polarization  P is  given 
by 


P - <N^  p12>  = N P21<v1v2*>  (1.9) 

where  N is  the  particle  density  and  the  brackets  indicate  an 
average  over  all  two  level  systems.  One  obtains  from  (1,1)  and 
(1.11)  by  dot-multiplication  of  both  sides  by  p12/itt 


_c_ 

Pl  2 ’ S' 

+ 9 

Pl  2 • 2 

ft2  8x 

a2at 

1 i*  . 

<v,  v 


12 


where 


2 fie 


(1.10) 


(1.11) 


The  equation  for  the  electric  field  (1.10)  completes  the 
system  of  equations;  the  solution  of  the  Zakharov-Shabat  equation 
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appears  directly  as  a drive  in  the  equation  of  the  field.  The 
system  of  equations  is  nonlinear  in  that  the  drive  is  nonlinear 

in  vi'  v2** 

Through  the  use  of  the  normalized  variables 


2c  = 6/fl 

z = ftx/c , 


the  equations  (1.6),  (1.7)  and  (1.10)  assume  the  form 


- - <2viV> 


3t 

9v, 

i 

3t 


+ iC  v.  = 


- t v. 


- ic  v. 


= - - (*  v. 


(1.12) 

(1.13) 

(1.14) 


These  equations  are  already  in  one  of  the  standard  forms  of 
inverse  scattering  theory.  McLaughlin  and  Corones  (1974)  have 
pointed  out  the  relation  between  the  linear  (Zakharov-Shabat ) 
problem  and  the  quantum  mechanical  equations  of  the  Josephson 
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junction.  They  also  touched  on  the  problem  of  S.I.T.  without 
making  the  connection  of  the  v's  with  the  wave  function  ampli- 
tudes. 

For  later  reference,  and  to  make  connection  with  the  Bloch 
equations,  we  also  list  the  differential  equations  for  the  density 
matrix  elements  (1.8).  They  follow  directly  from  (1.6)  and  (1.7) 


After  introduction  of  the  variable  t and  the  definitions 

X E 2pl2  = 2v1V2* 

N = P11  p22 

(1.17) 

(1.18) 


one  obtains  the  normalized  Bloch  equations: 

— + 2ia  = EN 
j 3t 

3JN 
3t 


= - - (E*X  + EX*). 
2 


3 
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Lamb  (1973)  used  the  Bloch  equations,  and  the  field  equation 

^.  = <X>  . (1.19) 

3z 

as  the  defining  equation  of  S.I.T.  Lamb  had  to  go  through  a series 
of  transformations  to  derive  the  Zakharov-Shabat  equations.  Our 
way  of  deriving  the  equations  shows,  much  more  simply,  that  the 
Bloch  equations  are  implied  by  the  Zakharov-Shabat  equations. 


T 
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II.  The  Scattering  Problem 

We  have  shown  that  the  nonlinear  self  induced  transparency 
equation  is  cast  naturally  in  terms  of  a set  of  linear  differ- 
ential equations  for  the  amplitudes  of  the  wave  functions  of  the 
upper  and  lower  levels  coupled  by  the  electric  field.  These  were 
the  equations  of  Zakharov  and  Shabat  (1972)  central  to  their  formu- 
lation of  the  inverse  scattering  problem  for  the  nonlinear  Schroe- 
dinger  equation  and  derived  by  Lamb  (1973)  from  the  density  matrix 
equations  of  the  two  level  system  by  a set  of  variable  transforma- 
tions. In  this  section  we  shall  elaborate  on  the  significance  of 
the  Zakharov-Shabat  equations. 

We  consider  them  to  be  a set  of  equations  of  inode  coup)  ing 
in  space , treating  T as  if  it  were  the  distance  coordinate,  <; 
the  propagation  constant  (c  is  real  by  definition);  the  amplitudes 
V1  and  v2  are  then  wave  amplitudes.  The  function  £.(t)  plays 
the  role  of  the  coupling  coefficient.  in  the  absence  of  an  £- 
field , 

v1  « exp  - id  (2.i) 

and 

v2  « exp  + id.  (2.2) 

1 

To  wav o Vj  propagates  in  the  -r  direction,  the  wave  v,  in 


the  +t  direction  (assuming  the  physicist's  definition  of  phase- 
delay  as  represented  by  the  factor  exp  i^x  (£  > 0)).  The  ampli- 
tudes v ^ and  are  functions  of  £ and  t.  Consider  briefly 

the  Fourier  transform 


d£  e i?y  v-^C,  t)  = V1(y,  x)  . 

< — oo 


(2.3) 


If  C is  taken  to  be  a propagation  constant,  £ = w/u  with  u 
the  phase  velocity  of  the  uncoupled  wave,  w the  frequency,  then 
y may  be  interpreted  as  a time  variable  (y  = ut) . This  further 
interpretation  endows  the  waves  and  with  dispersion-free 

propagation  at  group  velocity  u,  in  the  absence  of  t.  The 
original  self-induced  transparency  problem  involving  interactions 
of  electromagnetic  pulses  with  the  nonlinear  medium  requires  that 
£ (x)  has  to  vanish  at  |x|  ■*  °°.  Hence,  equations  (1.13)  and  (1.14) 
describe  coupling  of  waves  in  an  interaction  region  extending  from 
-oo  < r < +oo  , with  vanishing  interaction  in  the  limit  |x|  -*■ 

In  the  region  where  £.  ^ 0,  the  forward  and  backward  waves 
are  coupled.  The  coupling  is  lossless  in  the  sense  that  (for  real 

~ (|v1(C,  t ) | 2 + |v2(c,  t)  | 2 ) = 0 (2.4) 

> dx 


as  can  be  demonstrated  easily  from  Eqs.  (1.13)  and  (1.14).  In 
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(2.4)  waves  (1)  and  (2)  may  be  assigned  powers  Iv^tc,  t) | 2 and 
|v2(C»  x) |2.  According  to  (2.4),  both  waves  carry  power  in  the 
same  direction — say  the  +x  direction.  Because  they  have  oppo- 
sitely directed  group  velocities,  their  energies  must  be  of  oppo- 
site sign  (Pierce,  1974). 

The  concept  of  negative  small  signal  energy  is  widely  used 
in  plasma  physics  (Sturrock,  1961).  Negative  energy  commonly 
occurs  in  energy  conservation  principles  derived  from  the  linearized 
equations  of  motion  of  a nonlinear  system  which  contains  an  energy 
"reservoir"  (such  as  the  kinetic  energy  of  a moving  plasma  or  an 
electron  beam  (Pierce,  1974)).  Excitation  of  a wave  (usually  a 
socalled  slow  wave  (Sturrock,  1961))  may  lower  the  overall  energy 
of  the  system,  a fact  that  manifests  itself  in  terms  of  a negative 
energy  attributed  to  the  wave.  The  energy  is  quadratic  in  the 
excitation  amplitude  of  the  wave.  If  a negative  energy-wave  is 
coupled  to  a positive  energy  wave,  both  wave  amplitudes  may  grow. 

The  growth  of  positive  energy  is  balanced  by  the  growth  of  the 
negative  energy,  net  small  signal  energy  is  conserved.  One  example 
of  such  a system  is  the  Backward-Wave  oscillator  (Kleen,  1958). 

More  familiar  may  be  the  example  of  parametric  interaction 

(Yariv,  1976)  of  two  waves  of  frequencies  and  with  a 

pump  wave  of  frequency  w , so  that  u>,  + u>0  = w . The  phase 

P 1 2 p 1 

matching  condition  of  the  (col inear)  propagation  vectors  k^  and 
k 2 then  kp  = k^  + kj*  *n  s^-eac^y  state,  when  phase  matching 

is  not  realized,  one  may  define 
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C = — [ (k.  + k.)  - k ] 

y 1 2 P 


The  equations  of  parametric  coupling  between  the  two  waves  as  a 

function  of  the  spatial  coordinate  (t)  are  of  the  form  of 

* 

(1.13)  and  (1.14)  when  the  "fast"  spatial  variations  of  the  waves 
are  removed  and  only  the  slowly  varying  variations  of  "envelopes" 
are  considered.  The  energy  densities  must  be  reinterpreted  as 
photon  number  densities,  the  power  flows  as  photon  number  flows 
and  energy  conservation  as  photon  number  conservation.  To  be 
more  specific,  in  a parametric  process  of  the  type  where  a pump 
photon  produces  a "signal"  photon  at  frequency  u1  and  an 

"idler"  photon  at  frequency  u2 , the  number  of  signal  photons 
generated  either  spontaneously,  or  by  induced  emission,  must  be 
equal  to  the  number  of  idler  photons.  The  wave  amplitudes,  v^ 
and  v2  may  be  so  normalized  that  | v1  | 2 and  |v2|2  are  pro- 
portional to  the  number  of  photons  per  unit  length  in  the  inter- 
acting wave  (1)  and  (2).  Then  (with  t taken  as  the  distance 
variable) 


IvJ*  ± ± |v2|’  . 0 


are  the  Manley  Rowe  conservation  relations  applied  to  this  para- 
metric process  (Manley  and  Rowe,  1959;  Weiss,  1957).  A parametric 
instability  occurs  with  the  + sign  in  the  above  equation,  when 
the  group  velocities  of  waves  (1)  and  (2)  are  oppositely  directed. 

The  waves  v^  and  v2  have  further  properties  somewhat 


i 


I 


i 

I 
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analogous  to  lossless  coupling  of  electromagnetic  waves  that  obey 
reciprocity  relations.  Indeed,  from  (1.13)  and  (1.14)  it  is 
easily  shown  that,  given  a solution  (f  is  treated  as  a column 
matrix  of  components  f^  and  f2) 


f 


then 


( 1 


f 


* 

2 


‘V 


is  also  a solution  for  the  same  c (if  £ is  real).  Further, 
these  two  solutions  are  physically  different.  Indeed,  have  f 
describe  the  coupling  of  wave  v2  to  v^  via  £(x)  with  boundary 
conditions  as  indicated  schematically  in  Fig.  25a.  The  solution 
f is  the  one  shown  in  Fig.  25b.  The  function  f in  relation  to 
f is  like  the  time  reversed  solution  of  electromagnetic  waves 
used  to  demonstrate  reciprocity. 

One  may  generalize  (2.4)  to  show  conservation  of  "cross  power", 
i.e.  prove  the  conservation  law 


dx 


Uig 


lyl 


+ f2V> 


0 


(2.5) 


where  f and  g are  any  two  solutions  of  (1.13)  arid  (1.1!). 


Usi  ng 
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the  property  that 


is  a solution,  if  g is  one,  (2.5)  becomes 

— (fxg2  " f2^i)  = °*  (2-6) 

dx 

This  is  known  as  conservation  of  the  Wronskian  (Gardner  et  al.,  1967). 

Thus  far  we  have  studied  general  properties  of  the  scattering 
problem.  One  may  use  (1.13)  and  (1.14)  to  find  solutions  v^  and 
v2  for  given  £.  More  relevant  to  the  solution  of  the  self-induced 
transparency  problem  is  the  inverse  scattering  problem:  given  v^ 
and  v2 , what  coupling  function  F (t)  produces  this  particular 
v^  and  v^. 

The  S.I.T.  problem  calls  for  a very  special  kind  of  solu- 


tion v^ ( £,  t) , v2(S,  x).  Indeed,  £ is  the  parameter  describing 

the  detuning  of  the  two  level  systems  from  the  carrier  frequency 


wQ.  If  there  is  to  be  no  loss,  |v1(C»  t)  | = 1,  |v2U,  x)  | ® 0 

for  x ->  i.e.  every  two-level  system  has  to  start  from  the 

ground  state  .before  the  arrival  of  the  pulse  and  must  return  into 


the  ground  state  after  passage  of  the  puls^,  requirement  in 

„ ,hL,  „ * ■■ 

turn  calls  for  a scattering  , well  which  prpdUOtlB  hd  reflec- 


tion  (V  - 0)  for  an  in^ffienft  ’,fWav%'’  v,  (c,  x)  for' any  'c  1 

cnvu  oiTT*  1 bo-a^ac.  ot  W . »n«F 
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There  are  wells  that  are  capable  of  doing  this.  The  Schroedinger 
equation  for  a secant  hyperbolic  well  has  a continuum  of  eigen- 
states that  are  traveling  waves  outside  the  well  and  experience 
only  a phase  shift  as  they  pass  through  the  well  (Morse  and  Fesh- 
bach,  1953).  In  the  next  section  we  set  up  a method  for  obtaining 
the  shapes  of  reflection-free  wells  for  the  Zakharov-Shabat  equa- 
tions. 

% 
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III . The  Inverse  Scattering  Problem 


We  have  pointed  out  that  the  solution  of  self-induced  trans- 
parency calls  for  the  determination  of  a reflection-free  well.  The 
inverse  scattering  theory  determines  the  shape  of  a scattering  well 
from  the  scattering  data.  The  requirement  that  there  be  no  reflec- 
tion is  sufficient  to  find  shapes  of  reflection-free  wells. 

We  shall  develop  the  method  of  inverse  scattering  by  simple 
physical  reasoning.  Consider  a given  scattering  well  E(x)  with 
the  particular  solution  f(C,  t)  that  approaches  the  limit  as 


-f  ->  — oo 


lim  f(c,  t)  = 

T -►’-oo 


fil 

0 


-Ut 


(3.1) 


The  function  f(c,  t)  represents  an  experiment  in  which  a 
wave  v^  is  incident  from  x -*■  +°°,  partly  coupled  to  the  re- 
flected wave  v2  and  partly  transmitted. 

Next  consider  f defined  by 


which  is  obtained  by  time  reversal  of  f.  It  is  independent  of  f 


because  it  is  the  solution  to  a different  boundary  value  problem. 

Finally,  introduce  a third  solution  g(£,  t)  defined  by  its 
limit  at  x -*•  +°°: 


lim  g U , x) 

X -»-oo 


(3.2) 


Figure  26  shows  the  experiment  represented  by  g.  A mode  (wave) 
is  incident  from  the  left  upon  the  interaction  region,  partly  re- 
flected (coupled  to  the  backward  wave  v^)  and  partly  transmitted. 
g(C,  x)  may  be  written  as  a linear  superposition  of  the  two  inde- 
pendent solutions  of  the  (second  order)  Zakharov-Shabat  equation 
system. 

gU,  x)  = a ( £)  f(£,  x)  + b(U  ?U,  x).  (3.3) 

The  inverse  scattering  theory  enables  one  to  construct  the  scat- 
tering well  from  the  information  on  b(C)  and  a(C).  We  shall 
now  go  through  the  derivation  in  a "physical"  way.  Consider  the 
scattering  experiments  represented  by  g,  f,  and  f.  Remember 
£(<=£>  co/v)  is  interpreted  as  a propagation  constant  of  a disper- 
sion free  wave  and  hence  is  proportional  to  the  frequency.  g(£,  x) 
is,  at  first,  defined  on  the  real  ?-axis — for  a steady  state 
scattering  experiment.  For  a given  well,  the  function  gU,  x) 
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may  be  continued  analytically  into  the  upper  naif  of  the  complex 
^-plane.  This  can  be  interpreted  simply  in  terms  of  an  experiment 
using  a wave  source  (of  v2)  on  the  left  side  of  the  well  with  an 
amplitude  growing  exponentially  with  time.  Each  point  in  the  upper 
half  of  the  £ -plane  may  be  thought  to  have  associated  with  it  the 

space  (t-)  dependence  of  the  incident  transmitted  and  scattered 
waves.  We  show  such  plots  schematically  in  Fig.  27.  Fig.  27a  shows 
g(C,  t)  for  a point  in  the  upper  half  of  the  £-plane.  Shown  are 
the  space  (x-)  dependences  of  (a)  the  well  £(x),  (8)  the 

exciting  wave  g2  of  dependence  lim  exp  it;x,  and  (y)  the  reflected 

X4- 00 

wave  g which  is  caused  by  reflection  of  the  well  and  is  confined 
to  a finite  region  of  space. 

As  Im  r,  is  increased,  the  incident  wave  decays  more  steeply 

with  increasing  t,  and  as  |c|  -*•  “ the  reflected  wave  g^  van- 

ishes. Indeed,  the  faster  the  increase  with  time,  (the  greater 
Im  ;;)  , the  shorter  the  duration  of  the  interaction  and  the  weaker 
is  the  excitation  of  the  reflection  g In  the  limit  of  Im(c)  -*■  00 
the  reflected  wave  g^  decreases  to  zero.  The  same  can  be  said 
about  a spatially  very  rapidly  varying  wave,  |c|  very  large.  The 
reflection  vanishes  not  only  for  Im(C)  -*•  00  but  for  | i;  | -*•  °°  in 
general.  Fig.  5b  shows  an  analogous  schematic  sketch  of  the  exper- 
iment  described  by  f(£,  x)  a wave  incident  upon  the  well  from  the 

right.  Like  in  the  case  of  g(C,  x),  the  reflection  f2(S,  t) 

becomes  weaker  and  weaker  as  | £ | becomes  larger  and  larger  and 
vanishes  in  the  limit  | 1;  | 

l 


(r 
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If  the  well  F(t)  is  chosen  "deep  and  wide  enough",  the 
feedback  provided  by  the  well  will  permit  trapped  solutions  growing 
with  time  at  a characteristic  rate  Im  (>  0),  where  may 

assume  one  or  more  discrete  values,  depending  upon  the  depth  and 
width  of  the  well.  A trapped  solution  is  one  for  which  there  is 
no  external  excitation  g£  -*■  0 as  t -*  -°°.  The  excitation  within 
the  well  decays  away  from  the  well  in  both  directions.  With  a fixed 
amplitude  of  g ^ on  the  right  hand  side  of  the  well,  the  amplitude 


of  the  excitation  within  the  well  remains  finite  for  C -*• 

Indeed,  in  order  to  keep  the  excitation  amplitude  of  tne  well 
fixed  when  the  "resonance",  represented  by  the  trapped  solution, 
is  approached,  the  source  amplitude  has  to  decrease  until  it  vanishes 
at  C = C.  . The  trapped  solution  is  shown  schematically  in  Fig.  6. 
Power  is  carried  out  of  the  well  to  supply  the  increasing  energy 
storage  outside  the  well — the  energy  being  supplied  by  the  growth 
of  positive  energy  in  g2  at  the  expense  of  the  negative  energy  in 
g-^  via  the  instability  within  the  well.  The  same  applies  to 
f (£,  t)  and  indeed  g(£,  x)  and  f(C,  x)  cease  to  be  linearly 
independent  at  ^ because  they  both  describe  the  same  trapped 

solution . 

Finally  consider  F(C,  x)  . For  real  c, , the  function  is  the 
time  reversed  version  of  f (compare  Fig.  25b).  Before  we  con- 
tinue it  analytically  into  the  upper  half  plane  we  must  understand 
the  physical  situation  represented  by  the  solution  f.  According 
to  Fig.  25b,  waves  are  incident  from  both  sides  of  the  well.  The 
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well  is  excited  by  two  sources.  The  sources  are  so  phased  that  the 
"reflected"  wave  ?1  on  the  left  hand  side  is  suppressed.  This  is 
the  solution  that  is  to  be  continued  analytically  into  the  upper 
half  plane.  A sketch  is  shown  in  Fig.  27c.  In  the  lower  half 
plane  f(c,  x)  represents  an  excitation  of  the  well  by  two  sources 
decaying  exponentially  in  time,  again  phased  so  as  to  cancel  the 
wave  f^  as  t -*■  -°°.  (See  Fig.  27d.)  Trapped  solutions  exist  in 
the  lower  half  plane.  Indeed  solutions  of  the  Zakharov  Shabat 
equations  occur  in  complex  conjugate  pairs.  Therefore  f must 
contain  the  trapped  solution  decaying  in  time,  at  £ = £ * in  the 
lower  half  c-plane.  f2  is  constrained  on  the  left  hand  side  of 
the  well , when  £ approaches  z,y*  » no  singularity  occurs  in  f 
because  a decaying  trapped  solution  receives  power  from  the  energy 
in  the  collapsing  tails  outside  the  well  and  the  condition 

_ \ r it  t § 

f^*  -*•  e * in  the  limit  x -*•  -°°,  imparts  the  tail  with  a 
finite  amplitude. 

We  are  now  ready  to  exploit  the  scattering  (thought)  experi- 
ments represented  by  g,  f and  f,  through  the  use  of  complex 
function  theory.  The  function  g(C,  x)  exp  - i?x  is  well  behaved 
throughout  the  entire  upper  half  £ -plane  ; as  | C | ->  00  it  ap- 
proaches [^]  . b(z)  in  turn  approaches  lim  b(5)  = 1,  because 

_ Ub00 

in  this  limit  no  reflection  occurs  and  f becomes  equal  to  g 
asymptotically.  We  multiply  (3.3)  on  both  sides  by 


(exp  - icx)/bU)  . 


(3.4) 
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In  the  upper  half  plane,  the  left  hand  side  has  poles  at  C = ^ 
in  the  upper  half  plane  and  approaches  [^]  in  the  limit  | c,  | 
Because  the  left  hand  side  is  equal  to  the  right  hand  side,  the 
latter  has  the  same  poles  and  same  limiting  behavior  in  the  upper 
£ -plane. 

Next,  we  define  the  function  f(£,  x)  exp  - icx  in  the  lower 
half  plane.  This  function  in  turn  approaches  [^]  as  |c|  -*■  », 
and  has  no  singularities.  The  discontinuity  between  the  two 
functions  in  the  two  half  planes  on  the  real  C-axis,  £ = £ , 
is: 


a (£) 

f(5,  x)  exp  - i£x. 

b(C) 

We  appeal  next  to  complex  function  theory.  If  a complex  function 
is  well  behaved  throughout  the  entire  complex  £ -plane  with  the 
exception  of  a finite  number  of  poles  and  a discontinuity  on  the 
real  axis  then  the  function  can  be  written  in  terms  of  a sum  of 
contributions  of  the  poles  and  in  terms  of  an  integral  representing 
the  discontinuity  along  the  real  axis.  This  property  of  complex 
functions  is  in  one  to  one  correspondence  with  potential  theory 
in  two  dimensions.  If  a two  dimensional  electrostatic  potential 
has  a finite  number  of  point  (line)  charge  sources  and  has  certain 
discontinuities  along  a curve  (surface) , these  discontinuities 
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being  produced  by  surface  sources,  then  the  potential  may  be  found 
using  the  appropriate  Green's  function. 

The  formal  evaluation  is  presented  in  Zakharov-Shabat  (1972). 
Here  we  shall  confine  ourselves  to  the  special  case  of  a reflection 
free  well,  a(£)  -*■  0.  Then 

9—'  ■T--  exp  - id  = f(£,  t)  exp  - id  (3.5) 

b ( C) 

everywhere.  Strictly  speaking,  a (5)  cannot  be  set  equal  to  zero, 
it  can  only  be  made  to  approach  zero,  because  otherwise  the  two 
sides  of  (3.5)  cannot  be  balanced  when  | £ | °°  as  can  be  seen 

from  the  asymptotic  behaviors  of  g(C,  t),  f(C,  t)  and  f(C,  t) . 

The  only  singularities  left  are  the  zeros  of  b(£)  where 
gUk,  1 ) becomes  proportional  to  ck  f(Ck,  t).  Suppose  g(C,  t)/ 
b(d  has  simple  singularities  at  Ck,  so  that  it  may  be  written, 
taking  advantage  of  the  fact  that  g(£k,  t)  “ f(Ck#  t)  at  C = &k, 


— i-L  exp  - iCi 
b ( r ) 


ck  f|rv  T) 

c - w 


exp 


lCkT  + 


(3.6) 


Then,  from  (3.5) 

N f-'k  fUk/  T)  (0 

l exp  _ iC  T = f(£,  T)  exp  - id  - 

k = 1 5 - 5k  (1 
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If  we  set  C * we  obtain  2N  equations  for  the  2N 


functions 


fl^j/  t)  and  f2*  Uj*,  x) 


N Ck  fl(Ckf  x) 


exP  - iC.T  = f*  (t  * t)  exp  - U *x  (3.7) 


k-1  Cj*  - Ck  ' ‘2  -j  - • - -~j 


? W<V>  n 


k ■ 1 S'  v 


exp  + i?k*T  * 


tty  t)  exp  + Ujt  - i. 


(3.8) 


The  unity  in  Eq.  (3.8)  takes  care  ofthe  fact  that  exp  - i£x  t^tt,  x) 
approaches  unity  for  | £ | ■+<*>. 

From  f one  may  evaluate  the  scatter-free  well  ►.  ( t ) by 
comparing  results  in  the  limit  C - ••  Prom  (1.14)  and  the  assumed 

asymptotic  behavior  of  £ (Fig.  27b)  we  have  in  the  limit  of  large 
C = K (real) 


3f. 


3x 


- i5f,  = 


- - f * f.  - - i.  f*  e”iCx 
2 1 2 


or 


1 f°° 

f 2 exp  i£x  f*(t)  exp  - 2i£x  dx  = - _L_ 

2 4i£ 


exp  - i£x  £ * (x  ) 


Thus 


Mx)  = - 4i£  exp  - i£x  iim  f *(£ , r ) . 


" - 


H . & ) 
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From  (3.7)  we  evaluate  t)  in  the  limit  of  large  5 , by 

replacing  C by  £: 

N C f (C.  , x) 

f *(£,  x)  exp  - i£x  = l exp  - i£.x  (3.10) 

* k = 1 K K 


and  thus 


N 

Mx)  = - 4i  I Ck  f i T)  e^P  - i5Rx. 


(3.11) 


In  concluding  this  section  one  should  mention  one  subtlety 
of  the  approach  that  is  often  not  explicitly  brought  out,  but  which 
our  approach  helps  to  illuminate.  We  have  pointed  out  that  f(£,  x) 
exp  - Lex  and  f(Q,  x)  exp  - ii.x  do  not  behave  properly  in  the 
complex  plane,  but  only  in  concert,  in  the  sum  (3.5).  To  cancel  the 
singularities  of  the  individual  terms  the  coefficient  a/b  is  essen- 
tial. When  analyzing  the  reflection-free  well,  although  a(c)  -*■  0 
on  the  real  axis,  we  have  not  set  aU)  equal  to  zero,  because  that 
would  make  it  strictly  zero  in  the  entire  complex  plane  and  lead  to 

i 

a breakdown  of  the  analytic  continuation.  Physically,  this  means 

that  the  conditions  of  causal  excitation  as  utilized  in  construction 

— ■ - - - - ■ - - - - 

of  the  solutions  of  a scattering  well  and  a reflection-free  well 
are  compatible  only  in  a limit,  the  limit  a -*■  0. 
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IV.  Invariance  of  x,  with  z 

In  the  preceding  sections  we  have  reviewed  the  solution  of 
the  Z akharov -Shaba t equations  and  given  physical  interpretations 
to  the  mathematical  steps.  We  recognized  that  S.I.T.  implies 
the  existence  of  a reflection  free  "well"  G(x)  so  that  the  lower 
level  amplitude  v^(£,  x)  starts  with  unity  magnitude  at  x = +00 
and  returns  to  unity  magnitude  at  x -*■  -<*>.  Such  wells  exist  and 
one  of  them  is  the  well  known  secant  hyperbolic  to  be  reviewed  in 
section  VI.  It  is  not  obvious,  however,  that  such  reflection-free 
wells  fc(x)  are  consistent  with  the  propagation  equation  (1.12), 
which  has  not  been  used  as  yet.  When  Lamb  (1973)  solved  the  211- 
pulse  problem  using  the  Zakharov-Shabat  equations  he  did  not  prove 
consistency  in  general,  but  only  by  substitution  of  the  solution 
into  (1.12).  Later,  Ablowitz  et  al.  (1974a)  provided  a general 
proof.  The  proof  will  be  presented  here  in  a slightly  modified 
form. 

The  Zakharov-Shabat  equations  were  shown  to  be  the  direct 
consequence  of  the  two-level  systems  equations.  The  parameter  C 
was  the  detuning  parameter,  obviously  real.  When  solutions  were 
sought  for  "wells"  F(x)  of  given  reflection  characteristics  (E- 
fields  that  excite  the  upper  and  lower  level  by  a given  amplitude) 
the  c parameter  was  extended  analytically  into  the  complex  plane. 
This  meant  the  study  of  responses  of  the  well  to  "waves"  v^U,,  x) 
and  v2(C/  x)  excited  by  exponentially  growing  the  decaying  sources. 
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No  physical  connection  can  be  made  between  the  responses  v^tc,  t) 
and  v2(c,  x)  and  the  two  level  systems.  Indeed,  for  complex  £ 
the  responses  grew  exponentially  with  | x | . This  finds  no  counter- 
part in  the  two  level  system  for  which  the  normalization  holds 
| v^ | 2 + | v2  | 2 = 1,  a normalization  not  maintained  for  complex  5. 
This  means,  however,  that  the  analytic  continuation  into  the  Q- 
plane  of  the  Zakharov-Shabat  equations  is  a mathematical  procedure 
which  cannot  be  interpreted  physically  in  terms  of  the  two  level 
system,  even  though  the  parametric  oscillator  analog  permits 
"physical  pictures"  to  be  associated  with  the  mathematical  pro- 
cedure. The  complex  r,  has  a new  significance,  not  directly 
apparent  in  the  two  level  system  equations  or  the  Bloch  equations. 

A reinterpretation  of  these  equations  is,  therefore,  in  order. 

The  Zakharov-Shabat  equations  for  any  given  complex  £ define 
the  well  F (x)  in  terms  of  v^j;,  x)  and  v2(c,  x).  Conversely, 
a given  well  F(x)  and  boundary  conditions  at  infinity  define 
v1(^,  x)  and  v2U,  x) . The  propagation  of  the  well,  i.e.  its 
x,  z dependence  are  then  reflected  in  the  x,  z dependences  of 
v ^ and  v2.  Nonlinear  differential  equations  that  are  soluble 
by  the  inverse  scattering  method  have  the  remarkable  property 
in  common  that  the  solution  of  the  scattering  problem  for  a 
particular  value  of  one  of  the  independent  variables  (z  in  our 
case)  implies  the  solution  for  all  its  values.  Mathematically 
this  property  is  demonstrated  by  showing  that  c is  independent 
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I 

this  variable  (2).  Ablowitz  et  al.  proceed 

as  follows. 

1 

They  assume  that  v^^  and  v2  evolve  with  2 as  follows. 

1 

• | 

1 

3vl 

1 

“ A(£/  2 , t)v.  + B((,  z,  t)v~ 

3 z 2 

(4.1) 

i 

M 

: | 

i 

3v2 

: = T)V.  + DU,  2,  T)Vj. 

3 2 A *• 

(4.2) 

1 

"v 

i 

In  order  that  (4.1)  and  (4.2)  be  consistent  with 

(1.13)  and 

• 

II 

I II 

i 

(1.14),  cross  derivatives  31/3t3z  must  be  equal 

to  3 2/3z3t. 

ij 

ii 

Prom  the  requirement  that  C be  independent  of 

z one  finds 

• • [I] 

i 

■ 

A = -D 

(4.3) 

, I - 

where  an  integration  constant  has  been  set  equal 

to  zero,  and 
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These  equations  must  correspond  to  the  equations  of  the  system, 
if  the  hypothesis  of  invariance  of  c with  z is  to  be  proven 
correct.  To  cast  the  equations  of  the  system  into  the  form  of 
(4. 5) -(4. 7)  we  take  a Hilbert  transform  of  the  Bloch  equations 
(1.17)  and  (1.18)  in  the  following  manner.  The  parameter  c used 
in  those  equations  is  the  real  detuning  parameter.  To  distinguish 
it  from  the  generalized  £,  we  use  the  symbol  5 for  it.  Further, 
we  introduce  the  line  shape  function  of  the  two  level  systems  g(£) 
normalized  so  that 


g (£)  d£  = 1. 


The  Hilbert  transform  of  any  function  f(£)  with  respect  to  the 
line  shape  function  g(£)  is  then  defined  as 


H (f ) = i-  f 

4 J -o 


f(U  g(S) 
C - S 


where  £ is  assumed  to  be  in  the  upper  half  plane  and  proper 
indentation  of  the  contour  is  implied  when  ( approaches  the  real 
axis.  By  transforming  (1.17)  and  (1.18)  and  the  complex  conjugate 
of  (1.18)  we  find 


— = - (I  ?C  - - P*B 
3t  ^2  2 


(4.7) 


where  A,  B,  and  C are  the  Hilbert  transforms  of  N(£,  z,  x) , 
A(£,  z,  x)  and  \*  (£,  z,  C)  respectively.  Comparison  of  (4.7), 
(4.8)  and  (4.9)  with  (4.4),  (4.5)  and  (4.6)  shows  that  the  two 
sets  of  equations  can  be  made  identical  through  the  use  of  the 
field  equation: 
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V.  The  Spatial  Dependenc e 


In  the  preceding  section  we  have  shown  that  the  complex  para- 
meter £ is  an  invariant  (independent  of  z)  . The  solution  of 
the  problem  being  described  fully  by  the  constants  and  the 

eigenvalues  ck . The  latter  being  z -independent , the  z-dependence 
is  contained  entirely  in  the  constants  Ck* . The  C^'s  in  turn  are 
independent  of  x ! Using  this  fact,  equations  for  Ck(z)  maY  be 
obtained  very  simply  in  the  limit  x -*■  -°°.  We  need  to  evaluate 
f^C,  x)  and  f2*(C,  x)  and  then  use  the  field  equation  (1.12). 

In  this  limit,  we  find  from  (3.1)  for  f^(£,  x) 


lim  f1  (s,  x)  = e"1^1. 


(5.1) 


f2 t)  follows  directly  from  (3.5)  in  terms  of  f^c^,  x) 


N C, 


lim  f_*(£,  x)  exp  - i£x  = £ 


k ■ 1 5 - 


exp  - 2i£kx.  (5.2) 


Next  we  need  the  asymptotic  expression  for  £ from  (3.11) 


lim  £ = - 4i  £ C.  exp  - 2.i£, x. 
x->- °°  k * 1 K 


(5.3) 


100 


Introducing  (5.3),  and  vx  * fx  of  (5.1),  v2  = f2  of  (5.2), 
into  (1.12)  we  get  the  following  differential  equations  for  the 

Vs* 

N 3Ck  N i . 

E exp  - 2ickx  = Z — C,  < > exp  - 2U-t. 

k = 1 3z  k k»12*E-r  k 

sk 

(5.4) 


Equating  terms  of  the  same  t dependence  one  obtains  a simple 
differential  equation  for  CR  which  has  the  solution 


Ck(2> 


ck(°) 


(5.5) 


This  result,  along  with  (3.11),  allows  one  to  construct  the  £- 
field  for  any  assumed  distribution  of  poles.  We  consider  first 
the  general  case  and  then  apply  the  result  to  the  2ir-pulse  and 
zero-pi  pulse. 


Combining  (3.7)  and  (3.8)  we  obtain  an  equation  for 
exp  iCjT  which  is 


t) 


Z 

j 


fl«j' 


t)  exp  iCjT  « 1 


(5.6) 


with 


C *(Z)  C.(z) 

D..  = 6..  + l exp  i(2Cv*  - C-  - £.)t.  (5.7) 

13  13  k (C-  - Ck*)  (CR*  - Cj)  K i D 


When  (5.6)  is  inverted  and  introduced  into  (3.11)  we  find 


F ( t , z ) = - 4 i Trace [MD  x ] 


(5.8) 


where 


-i(CH  + ?,)t 
Mij  = C j (z)  e 3 . 


(5.9) 


Equation  (5.8)  looks  like  Eq.  (16)  of  Lamb  (1973)  but  his  defini 
tions  of  M and  D are  wrong. 
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VI.  Solution  for  2ir-Pulse 

The  number  of  assumed  singularities  determines  the  complexity 
of  the  solution.  We  start  by  assuming  the  existence  of  only  one 
pole  j and  find  the  corresponding  solution  from  (5.8) : 


F(x)  = 4iC . g-xP  ~ 2(iotT  ~ 

3 |cJ! 

1 + — J — exp  4gt 

|23|2 


(6.1) 


From  (5.5)  we  find  the  z-dependence  of  C.  to  be 


C (z)  = C . (0)  exp  i > . 


2 E,  - a - i8 


(6.2) 


Introduction  of  this  expression  into  (6.1)  gives 


£(t,  z)  = 48  exp  - 2ia(x  - x^)  sech  28(x  - xq) 


where 


ax,  = — < £ a > 

4 ( £ - a)  2 + 82 


(6.3) 


0xo  = - < 


4 ( c “ a) 2 + B2 


(6.4) 
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C,  (z)  exp  - 2iC.  x =*  CK(0)  exp  - 2i(-a  + i3)x  exp  — < ^ — St) — i_L§_>  . 
b b b 2 U + a)2  + 32 

(7.2) 

The  speeds  of  the  functions  responsible  for  the  envelopes  are 
obtained  from 

23dx  = — < - > dz 

2 ( £ + a)  2 + 3 2 

respectively. 

The  speeds  will  be  different,  unless 

< 1 > = < 1 > (7.3) 

( £ - a)  2 + 32  ( £ + a)  2 + 32 

which  is  true  only  for  a symmetric  line  function  centered  with  the 
carrier  frequency  of  the  pulse.  When  the  line  is  symmetric,  then 
the  speed  of  "phase  propagation"  determined  from 

1 

± 2adx  = - < L_L_“ > dz 

2 ( E,  + a)  2 + 32 

is  also  equal  for  both  constituents.  In  this  case,  one  finds  from 
(5.8)'  after  some  algebra 


{ 
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Mt,  z) 


6 


= 88 


cos  2a(x  - x. ) cosh  2B(t  - x ) - — sin  2a(x  - x, ) sinh  28(x  - x ) 
l ° a l 1 

O 2 

cosh223(x  - x ) + — sinz2a(x  - x, ) sinh22a(x  - x, ) 

° a2  1 1 


(7.4) 


where 


Bx. 


ax, 


i z <- 


U - a)  2 + 8: 


l z <- 


S - a 


U - a)2  + B: 


(7.5) 


(7.6) 


This  is  the  same  result  in  a simpler  form  than  the  one  obtained 

by  Lamb  (1974 ) . 

when  the  poles  are  not  picked  symmetrically  with  respect  to 
the  ImU)  axis,  new  solutions  may  be  generated.  In  particular, 
one  is  interested  in  finding  solutions  that  do  not  split  up  into 
2 tt  pulses  and  possess  a spectrum  not  centered  with  the  material 
line.  Whether  such  solutions  exist  may  be  decided  with  the  aid 
of  (5.8).  In  limit  x -*■  °°,  the  expression  for  E(z,  x)  contains 
x - z dependences  only  in  the  combination  C^(z)  exp  - 2i£^x  or 
its  complex  conjugate.  This  combination  gives  a very  special 
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propagation  velocity  of  the  envelope  exp  20^t.  The  propagation 
velocity  (or  rather  its  inverse)  is  given  by  (compare  7.5) 


At 

Az 


1 <- 


a - cu)2  + e.- 


->  . 


If  the  solution  is  to  describe  a pulse  that  does  not  break  up, 
the  envelope  speeds  of  C exp  20  x and  C.  exp  20.  x must  be 

a a O D 

the  same.  This  requires 


< - > = < - > . (7.7) 

U - «a)a  + 6/  <e  - «b)2  + 0b2 

If  we  are  considering  a sharp  line,  with  detuning  £ , the  above 
may  be  written 

“a’  + - 2«oaa  * V + 6b’  - 2?o<V  <7'81 

Equation  (7.8)  puts  the  eigenvalues  and  £b  on  a circle 

centered  at  £q.  The  exponentials  exp  - ic^x  obtained  for  the 
case  of  no  detuning,  £q  = 0,  are  simply  multiplied  by  exp  - i£Qx 
Equations  (5.7)  and  (5.8)  show  clearly  the  effect  of  such  a multi- 
plication; does  not  change,  is  multiplied  by  exp  - 

2if.Qi  — a simple  frequency  shift.  The  sped  rum  is  shifted  by 
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? 


precisely  the  right  amount  to  center  it  again  with  the  line  of  the 
material  medium.  The  assumed  detuning  has  been  removed  by  the 
requirement  that  the  two  envelopes  propagate  at  the  same  velocity. 


. rt<,  v r 
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Conclusions 

We  have  shown  that  the  Zakharov-Shabat  (Z-S)  equations  are 
identical  with  the  two-level  equations  of  the  S.I.T.  problem. 

In  order  to  be  aided  by  physical  intuition  in  the  general  solution 
of  the  equations  we  showed  that  they  are  analogous  to  the  spatial 
differential  equations  of  the  parametric  amplifier-oscillator 
problem,  where  the  ^-parameter  of  the  Z-S  equations  plays  the 
role  of  frequency.  The  "time"  variable  t of  S.I.T.  is  inter- 
preted as  a spatial  variable  in  the  analog  problem. 

The  analytic  extension  of  the  scattering  functions  into  the 
upper  and  lower  half  planes  and  their  properties  were  made  under- 
standable by  the  physical  intuition  developed  for  the  excitation 
of  a given  "twoport"  from  transmission  lines  connected  to  the 
twoport  by  exponentially  growing  or  decaying  sources.  The  occur- 
rence of  zeros  of  the  reflection  coefficient  in  the  upper  half 
C-plane  corresponded  to  the  growing  oscillations  of  the  parametric 
oscillator. 

The  analytic  continuation  of  the  scattering  equation  con- 
taining the  zeros  of  the  reflection  coefficients  as  poles  enabled 
one  to  find  the  shape  of  the  scatter-free  wells,  the  x-dependence 
of  the  E-field  of  S.I.T.  at  one  value  of  the  spatial  coordinate 
z of  the  original  S.I.T.  problem.  It  was  shown  (following 
Ablowitz  et  al.)  that  the  position  of  the  poles  is  z-independent . 
This  fact,  and  the  fact  that  the  residues  are  independent  of  the 
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time  variable  t was  used  to  get  simply  the  spatial  (z-)  de- 
pendence of  the  E-field.  2ir-  and  0-ir  pulse  solutions  were 
obtained  very  directly.  We  also  showed  that  there  is  no  soliton 
solution  corresponding  to  an  off -resonant  zero-pi  pulse. 


i 


i 
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Figure  Captions 


Figure  1 The  linear  absorption  over  1 km  path  in  the  atmosphere 
(due  to  I^O)  as  a function  of  frequency  (cm-^) . 

The  curves  were  calculated  by  S.  A.  Clough  using  AFGL 
Line  Parameter  Compilation. 


Figure  2 Flow  Chart  for  Inverse  Scattering  Method  applied  to 
S.I.T. 


Figure  3 


Figure  4 


Zero-pi  pulse  propagation  for  t = tf  = 100  psec  as 
a function  of  time,  with  space  as  parameter. 

(a)  normalized  E-field 

(b)  polarization 

(c)  population 

(d)  intensity  (watts/cm2) 

(e)  magnitude  squared  of  Fourier  spectrum. 

Three-dimensional  displays  of  zero-pi  pulse  propagation 
for  t = tf  = 10  psec  as  function  of  time  and  space. 

(a)  normalized  E-field 

(b)  polarization 

(c)  population 

(d)  intensity  (watts/cm2). 
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Figure  5 


Figure  6 

I 
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Figure  7 


Lossy  propagation  of  zero-pi  pulse  versus  distance 
(cm);  t = tf  = 100  psec;  values  given  in  (7.10). 

(a)  pulse  area 

l 9 

(b)  normalized  pulse  energy 

(c)  linear  and  nonlinear  absorption  coefficients 

(d)  ratio  of  nonlinear  absorption  coefficient  to 
linear  value. 

Lossy  propagation  of  zero-pi  pulse,  with  distance 
(cm)  as  parameter;  t = tf  = 100  psec;  values 
given  in  (7.10) . 

(a)  normalized  electric  field 

(b)  population 

(c)  intensity  (w/cm2 ) 

(d)  squared  magnitude  of  Fourier  transform. 

Lossy  propagation  of  zero-pi  pulse  versus  distance 
(cm);  t0  = tf  = 10  psec;  values  given  in  (7.12). 

(a)  pulse  area 

(b)  normalized  pulse  energy 

(c)  linear  and  nonlinear  absorption  coefficients 

(d)  ratio  of  nonlinear  and  linear  absorption  coef 
f icients . 
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Figure  8 Lossy  propagation  of  zero-pi  pulse  with  distance 

(cm)  as  parameter;  t = tf  = 10  psec;  values  given 
in  (7.12)  . 

(a)  normalized  electric  field 

(b)  population 

(c)  power  (w/cm2 ) 

(d)  squared  magnitude  of  Fourier  transform. 

Figure  9 Lossy  propagation  of  zero-pi  pulse  versus  distance 
(cm);  t = tf  * 1 nsec;  values  given  in  (7.14). 

(a)  pulse  area 

(b)  normalized  pulse  energy 

(c)  linear  and  nonlinear  absorption  coefficients 

(d)  ratio  of  nonlinear  and  linear  absorption 
coefficients . 

Figure  10  Lossy  on-  and  off-resonant  propagation;  tfi  = tf  = 

100  psec;  Fourier  transform  magnitude  squared  as  a 
function  of  frequency  (cm  1)  and  space  (cm) . 


(a) 

on-resonance  case 

(b) 

Aco  = 0.05 

cm 

(c) 

Aw  = 0.10 

cm  ^ 

(d) 

o 

CN 

• 

o 

. u 

3 

< 

cm  ^ 
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Figure  11  Three-dimensional  displays  of  propagation  for  off- 

resonant  pulse  versus  time  and  distance  t = t-  = 

e f 

100  psec;  Aui  = 0.05  cm  ^ 

(a)  normalized  electric  field 

(b)  polarization 

(c)  normalized  population 

(d)  pulse  intensity  (w/cm2 ) . 

Figure  12  Three-dimensional  displays  of  normalized  electric 

field  versus  time  and  distance;  t = t_  = 100  psec. 

6 t 


(a) 

A u 

= 0 

(b) 

Au> 

= 0.10 

cm  *■ 

(c) 

Aw 

= 0.20 

cm  * 

(d) 

Aui 

= 0.5  i 

cm 

Figure  13  Three-dimensional  displays  of  normalized  real  polari- 
zations (imaginary  parts  are  not  shown)  versus  time 
and  distance;  t = tf  = 100  psec. 


(a) 

Au>  = 

0 

(b) 

Aoi  = 

0.10 

cm  * 

(c) 

Aui  = 

0.20 

-1 

cm 

Figure  14  Three-dimensional  displays  of  population  versus  time 
and  distance;  t = tf  = 100  psec. 
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Figure  15  Three-dimensional  displays  of  pulse  power  versus  time 


and 

distance; 

fce  " 

(a) 

Aw  = 0 

(b) 

Au>  = 0.10 

cm  ^ 

(c) 

Aw  = 0.20 

-1 

cm 

Figure  16  Lossy  propagation  versus  distance  of  on-resonant 
zero-pi  pulse  for  t = tf  = 100  psec. 

(a)  pulse  area  (real  A and  imaginary  B) 

(b)  normalized  pulse  energy 

(c)  deviation  of  the  pulse  velocity  u from  the 
speed  of  light  in  terms  of  1 - u/c 

(d)  loss  coefficients  for  linear  and  nonlinear 
absorption 

(e)  the  ratio  of  nonlinear  to  linear  loss. 


Figure  17  Lossy  propagation  versus  distance  of  off-resonant 

(Ago  = 0.05  cm  *)  zero-pi  pulse  for  t = tf  = 100  psec. 

(a)  pulse  area  (real  A and  imaginary  B) 

(b)  normalized  pulse  energy 

(c)  deviation  of  the  pulse  velocity  u from  the 
speed  of  light  in  terms  of  1 - u/c 

(d)  loss  coefficients  for  linear  and  nonlinear 
absorption 

(e)  the  ratio  of  nonlinear  to  linear  loss. 


Figure  18  Lossy  propagation  versus  distance  of  off-resonant 

(Aw  = 0.10  cm  zero-pi  pulse  for  t = tf  = 100  psec. 

(a)  pulse  area  (real  A and  imaginary  B) 

(b)  normalized  pulse  energy 

(c)  deviation  of  the  pulse  velocity  u from  the 
speed  of  light  in  terms  of  1 - u/c 

(d)  loss  coefficients  for  linear  and  nonlinear 
absorption 

(e)  the  ratio  of  nonlinear  to  linear  loss. 

Figure  19  Lossy  propagation  versus  distance  of  off-resonant 

( Aw  = 0.15  cm  *)  zero-pi  pulse  for  t = tf  = 100  psec. 

(a)  pulse  area  (real  A and  imaginary  B) 

(b)  normalized  pulse  energy 

(c)  deviation  of  the  pulse  velocity  u from  the 
speed  of  light  in  terms  of  1 - u/c 

(d)  loss  coefficients  for  linear  and  nonlinear 
absorption 

(e)  the  ratio  of  nonlinear  to  linear  loss. 

Figure  20  Lossy  propagation  of  zero-pi  pulse  versus  distance; 

18.577  cm  ^ line;  t = tf  = 100  psec. 

(a)  pulse  area 

(b)  normalized  pulse  energy 

(c)  linear  and  nonlinear  absorption  coefficients 

(d)  ratio  of  nonlinear  to  linear  absorption  coefficients. 
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Figure  21  Lossy  propagation  of  zero-pi  pulse  versus  distance; 

20.204  cm  ^ line;  te  = tf  = 100  psec. 

(a)  pulse  area;  A,  B,  C,  D corresponds  to  |m|  = 

3,  2,  1,  0 

(b)  normalized  pulse  energy 

(c)  linear  (B,  D,  F,  H)  and  nonlinear  (A,  C,  E,  G) 
absorption  coefficients  for  |M|  = 3,  2,  1,  0 

(d)  ratio  of  nonlinear  to  linear  absorption  coef- 
ficients; A,  B,  C,  D corresponds  to  | M { =3, 

2 » 1,  0 

(e)  linear  and  nonlinear  absorption  coefficients 
averaged  over  all  M-levels 

(f)  ratio  of  averaged  nonlinear  to  linear  absorption 
coefficients . 

Figure  22  Lossy  propagation  of  zero-pi  pulse  versus  distance 

25.080  cm  ^ line;  tg  = tf  = 100  psec. 

(a)  pulse  area;  A,  B corresponds  to  |m|  =2,  1 

(b)  normalized  pulse  energy 

(c)  linear  (A,  B)  and  nonlinear  (c,  D)  absorption 
coefficients 

(d)  ratio  of  nonlinear  to  linear  absorption  coef- 
ficients A,  B corresponds  to  |M|  = 2,  1 


120 


Figure  23 


Figure  24 


Figure  25 


(e)  linear  and  nonlinear  absorption  coefficients 
averaged  over  all  M-levels 

(f)  ratio  of  averaged  nonlinear  to  linear  absorption 
coefficients . 

Lossy  propagation  of  zero-pi  pulse  versus  distance; 

30.560  cm  ^ line;  = tf  = 100  psec. 

(a)  pulse  area;  A,  B,  C corresponds  to  | M | = 2,  1, 

(b)  normalized  pulse  energy 

(c)  linear  (B,  D,  F)  and  nonlinear  (A,  C,  E) 
absorption  coefficients  for  | M | =2,  1,  0 

(d)  ratio  of  nonlinear  to  linear  absorption  coeffi- 
cients, A,  B,  C corresponds  to  |M|  =2,  1,  0 

(e)  linear  and  nonlinear  absorption  coefficients 
averaged  over  all  levels 

(f)  ratio  of  averaged  nonlinear  to  linear  absorption 
coefficients. 

Propagation  of  CH^F  line  in  presence  of  absorption 

of  18.577  cm  ^ line. 

(a)  energy  as  a function  of  distance 

(b)  the  absorption  coefficient. 

"Experiments"  defining  f and  f. 
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Figure  1.  Calculated  using  AFGL  line  parameters  compilation 

by  S.  A.  Clough 
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Figure  4d. 
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Figure  6a.  Normalized  electric  field  in  lossy  case 
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Figure  lib.  Lossy  Ott  pulse  propagation  Ao)  = 
0.05  cm  1j  normalized  polarization 
100  psec  = t = t , . 
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Figure  13b.  Lossy  Ott  pulse  propagation  t = 

t~  = 100  psec;  normalized  polarization 
' r -1 
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Figure  14a.  Lossy  Ott  pulse  pri 
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Figure  15c.  Lossy  Ott  pulse  propagation  t 

tf  = 10  0 psec;  pulse  power  Au> 

0.20  cm"^ . 
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. Lossy  propagation  of  zero-pi  pulse  versus  distance 
1R.C'7'’  cn  1 line;  t * tf  = 100  psec . 
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Figure  22.  Lossy  propagation  of  zero-pi  pulse  versus  distance 
25.080  cm  line;  t^  = t^  = 100  psec. 
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Figure  23.  Lossy  propagation  of  zero-pi  pulse  versus  distance; 
30.560  cm  line;  t « t^  = 100  psec. 
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(b)  normalized  pulse  energy 
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Figure  27  The  functions  g,  f and  f in  complex  t,  plane 
for  assumed  well  £(t).  Arrows  indicate  direction 
of  group  velocity.  The  reflected  wave  stops  abruptly 
at  end  of  well  as  indicated. 


